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Abstract 



In this work, the complete one loop calculation of meson-meson scattering amplitudes within U (3)(8> 
t/(3) chiral perturbation theory with explicit resonance states is carried out for the first time. Partial 
waves are unitarized from the perturbative calculation employing a non-perturbative approach based 
on the N/D method. Once experimental data are reproduced in a satisfactory way we then study the 
resonance properties, such as the pole positions, corresponding residues and their Nc behaviors. The 
resulting Nc dependence is the first one in the literature that takes into account the fact that the rji 
becomes the ninth Goldstone boson in the chiral limit for large Nc- Within this scheme the vector 
resonances studied, p{770), if* (892) and </)(1020), follow an Nc trajectory in agreement with their 

■ standard qq interpretation. The scalars /o(1370), ao(1450) and if*(1430) also have for large Nc a qq 
I pole position trajectory and all of them tend to a bare octet of scalar resonances around 1.4 GeV. The 

t:;^ ■ /o(980) tends asymptotically to the bare pole position of a singlet scalar resonance around 1 GeV. The 

OO ! cr, K and ao(980) scalar resonances have a very different Nc behavior. The case of the a resonance is 

■ analyzed with special detail. 

o 

^ ■ 1 Introduction 

J> A completely reliable theory to describe the physics in the intermediate energy region (resonance region) 
^ of QCD is still missing nowadays, since neither the low energy effective theory of QCD, namely chiral 
;h ■ perturbation theory (xPT) nor perturbative QCD is valid in this region. The appearance of very broad 
resonances, such as a and k [Ij, makes the discussion even more difficult and thus poses a great challenge 
for theorists. Several attempts have been done to address the question by implementing non-perturbative 
methods inspired from the S-matrix theory [2l[3l[l], typically on the xPT results. One has approaches like 
the Inverse Amplitude Method (lAM) [5l[6l[71[8] to establish the unitarized amplitudes, others are based 
on the N/D equations [9l [TTl [121 E], Roy equations [HI [151 IS]) or on building specific dispersion 
relations for the scattering amplitudes [171 118j . etc. Although the different approaches determine values 
for the pole positions of the broad resonances a and k in good agreement between each other, their 
nature is still a controversial subject. One way to get insight into them is to track the Nc trajectories of 
resonance poles in the complex plane. All the previous works along this research line [19 \ [20 l [2T1 \22 \ [23l I24j 
are performed within SU{3) or SU{2) |25j . In the former case the degrees of freedom correspond 

to the lightest octet of pseudoscalar mesons, pions tt, kaons K and the isoscalar r^g, while the singlet 
pseudoscalar rji is considered as a heavy field buried in the chiral counterterms. In the case of SU{2) 
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xPT only pions are taken as degrees of freedom. In Ref. |26j the scalar-isoscalar states, including the cr, 
are investigated within the large- A'^c Regge approach. 

C/(3)(8'C/(3) chiral symmetry in QCD is broken because of quantum effects that violate the conservation 
of the singlet axial vector current by the f/A(l) anomaly [27 ^ 128 ^ [29], As a result the singlet pseudoscalar 
T/i is not a pseudo-Goldstone boson [30j. Nevertheless, from the large Nc QCD point of view, the quark 
loop responsible for the Ua{^) anomaly [27] is I/Nq suppressed, thus indicating that the rji becomes 
the ninth pseudo-Goldstone boson in the large Nq limit |3H I32j. Hence it is necessary to include the 
r]i meson as a dynamical degree of freedom if one attempts to discuss the Nc trajectories of various 
resonances. This important fact was lacking in previous studies and is one of the main motivations of 
our current work. 

Thanks to the large Nq QCD, the singlet 7?i field can be conveniently incorporated into the effective 
field theory by enlarging the number of degrees of freedom of the theory from the octet of pseudo- 
Goldstone bosons to the nonet, which is usually called C/(3) xPT [33] . One advantage of foremost 
importance in our work is to use C/(3) for considering the running with Nc of resonance poles 

because then one has a framework that conceptually admits to consider the large Nc limit in xPT. E.g. 
it has then the proper number of degrees of freedom. 

Compared with SU{3) xPT, employed in the previous studies [El [20l EH [22l [23l [24] , a novel ingredient 
in U{3) xPT is the r/i mass term from the Ua{^) anomaly, which has nothing to do with the current 
quark masses. The appearance of this new scale in xPT could totally breakdown the well celebrated 
chiral power counting. The introduction of the 1/Nc expansion in U{3) xPT fixes the problem, since 
the singlet mass squared Mq behaves as 0{1/Nc) in the large Nc limit. Thus, in U{3) xPT there are 
three expansion parameters: momentum, quark masses and 1/Nc, giving rise to a joint triple expansion 
5 ^ ruq ~ 1/Nc- However, one should bear in mind that the value of the singlet mass Mq, mainly 
determined from the masses of the physical states rj and r]' , is not a small quantity (for a review see 
Ref. j34j ) . In connection with this problem Ref.[35j used the infrared regularization method (IR), which 
is proposed to cure chiral violating terms from loop corrections in baryon chiral perturbation theory in the 
low energy region \36\ 137] . to investigate the r]' physics within U{3) xPT. In IR the basic one loop scalar 
integral is divided into the infrared singular and regular parts. The regular part is not considered in IR 
because its effects are absorbed into the low energy constants of the theory. One only needs to consider 
the infrared singular part. However, contrary to baryon chiral perturbation theory, nothing prevents the 
appearance of large r]' masses in the vertices from derivatives acting on the external ij' fields, while this is 
not an issue in baryon chiral perturbation theory because of baryon number conservation. By the same 
token, in meson the total energy squared in the center of mass frame (CM) is not restricted to be 
around one, two, etc r]' masses, which could indeed spoil the validity of a loop calculation within IR. 
Although there are some special processes, such as rj' — )■ r/vrvr [35], which happens to be covered by the 
constrained energy region of IR, in general one can not naively apply IR in the pure meson sector. 

In the present discussion, we employ standard dimensional regularization, as in conventional SU{3) 
xPT [25], already applied in U{3) xPT to calculate the pseudoscalar weak decay constants in Ref. [38j . 
Then, the triple chiral expansion scheme is necessary for preserving the power counting. Investigations 
within U{3) xPT along this line mainly focused on the construction of higher order Lagrangians [39\ HQ] 
and the tree level and one- loop calculations for rj — r]' mixing |4H [38] . An important contribution of our 
work is to ofi^er the first complete one-loop calculation in U{3) xPT for meson-meson scattering, including 
both the loop diagrams contributing to the pseudo-Goldstone masses and decay constants, as well as the 
genuine ones of meson- meson scattering. A preliminary version of this calculation is given in Ref. [42]. 

Previous analyses from other groups only considered tree level amplitudes or performed partial one- 
loop calculations [43l HU [45] . The rj' rjTnr decay was studied in detail. Ref. [35j performed a one- loop 
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calculation within IR, and Ref. [1^ undertook a two-loop calculation within the framework of non- 
relativistic field theory. However, in this decay the three-momenta of the pions are not small compared 
with their mass for some region of the kinematics. Recently, Ref. [ITj discussed the same process within the 
triple expansion scheme by considering tree level amplitudes and part of the s-channel loops, resumming 
vrvr final state interactions. 

Since one of the main interests of our work is to study resonance properties, we will explicitly include 
bare resonance fields in the Lagrangian within the framework of resonance chiral theory [481 149] . However, 
including resonances as dynamical degrees of freedom is not enough to guarantee that one can safely apply 
the perturbative results up to the resonance region, since the loops of pseudo-Goldstone bosons, specially 
the unitarity or the s-channel loops, start to play an important role around the energy region where 
resonances emerge. Thus, a proper way to resum unitarity loops is crucial to study resonance properties. 
We use in this work the method provided in Ref. [10] to accomplish this resummation. Notice also that 
not all the resonances that stem in our study correspond to bare fields because new resonances come out 
when strong enough interactions are resummed. 

The present paper is organized as follows. Section [2] is devoted to the introduction of the relevant 
chiral Lagrangian. The structure of the perturbative scattering amplitudes is elaborated in Section [3l 
which is followed by the partial wave projection and its unitarization in Section |4l The phenomenological 
discussion, including the fit quality and features of resonances, such as masses, widths, residues and Nc 
behavior, is given in detail in Section [5l We conclude in Section [6l 



2 Relevant Chiral Lagrangian 

The chiral Lagrangian at leading order in U{2>) xPT reads 

= — K^'') + + — Mo'ln^detn, (1) 

where (. . .) denotes the trace in fiavor space and the last term corresponds to the Ua{^) anomaly rji mass 
term. The definitions for the chiral building blocks are 

X+ = u+xu^ + UX^U , 

rr 2 

Df,U = d^U-ir,,U + iUlf,, 

X = 2B{s + ip), (2) 

where ,1^ , s ,p stand for external sources and the pseudo-Goldstone bosons are collected in the matrix 



TT 



In addition, F is the axial decay constant of the pseudo-Goldstone bosons in the simultaneous chiral 
and large Nc limit. In the same limit the parameter B is related to the quark condensate through 
(0|^g-^|0) = —F^B5''K The explicit chiral symmetry breaking is implemented by taking the vacuum 
expectation values of the scalar external source field s = diag(mM, m^, m^), with mq the light quark 
masses. Throughout we always work in the isospin limit, i.e. taking rriu = m^. 
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We follow the framework in Ref . [48j to include bare resonance fields, where the interaction terms 
of pseudo-Goldstone bosons and resonances are invariant under chiral symmetry [50] and the discrete 
symmetries of charge conjugation C and parity P. Due to the presence of heavy resonance states, the 
momentum expansion is not valid any more in this theory. Nevertheless the 1/Nc expansion provides 
another guided principle to construct the Lagrangian. The generic Nq leading structure of the interacting 
operator in resonance chiral theory has only one flavor trace and in a schematic way can be written as 

~ {RiR2...RjX^^^), (4) 

where x^"^ denotes the chiral building block with chiral order 0{p'^), that only incorporates the pseudo- 
Goldstone bosons and the external source fields, and Ri stands for the resonance fields. Although in the 
resonance chiral theory the interacting terms with higher chiral orders (n > 4) are not suppressed in 
general, many of them are absent if one invokes the short distance constraints from QCD in large Nq |49] 
and does not allow fine tuning between different operators [5Tj. To be practical in the phenomenological 
discussion we restrict ourselves by including the full set of operators with n < 2. Concerning the number 
of resonance states in the interacting vertices, the most relevant ones in meson-meson scattering consist 
of one single resonance field. 

In regards to the operators attached only to the interactions between the pseudo-Goldstone bosons 
beyond leading order, it is commonly believed that they encode the high energy dynamics of the under- 
lying theory, which could be represented by the heavier states, such as resonances. It has been shown 
that at 0{p'^) in SU{3) xPT the low energy constants (LECs) are saturated in good approximation by 
the lowest multiplet of resonances and thus no additional pieces of LECs seem to be typically needed in 
the theory [48lli9]. 

In the present work, we exploit this assumption on the saturation of the LECs by the lightest reso- 
nances, so that, instead of local chiral terms contributing to meson- meson scattering we take the tree 
level exchanges of the scalar and vector resonances. In this way we keep all local contributions to meson- 
meson scattering up to and including 0{6^), while also generating higher order ones. We also calculate 
in addition the one-loop contributions that count one order higher in 6. 

The relevant resonance operators for the interactions with the pseudo-Goldstone bosons read |48] 



= CdiSsU^u'') + Cm{S8X+) + CdSl{u^U^) + CmSl{x+) 

+Cd{SQUf,){u^) + CmSiln'^ detu. 



(5) 
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where the resonance states are collected in the building block R = S^V , 
Si 
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The corresponding kinetic terms for resonance states read ^48j 



= l{V^^SsV,Ss-MlSl) + ^{d>^S,d,S,-MlSf), (9) 



where 



V^R = df,R+[r^,R], R = V,S, 

1 



(10) 



Compared with Ref. [IH] two additional l/Nc suppressed operators appear in Cs due to the inclusion 
of the singlet rji. These are the monomials proportional to q and Cm in Eq. ([5]). From the exchange of 
the scalar resonances these terms give rise to tree-level meson-meson contributions that are higher order, 
at least 0{5'^). In addition, these new operators mainly contribute to processes involving the r]' mesoru 
and since we deal with experimental data related to vr, K and rj in the present discussion, as shown 
explicitly below, the states with r]' only enter through an indirect way. So their effects are rather tiny in 
the current discussion and we discard these two new terms throughout. We have checked that if included 
our results barely change. For the remaining parameters of scalar resonances, instead of imposing the 
large Nc relations to the couplings and masses of the octet and singlet, such as Cd^m = Cd.ml^f^ and 
Mq-^ = Mo-g [Ml, we free them in our discussion. In this way, we consider effects beyond leading order of 
'i/Nc in the scalar resonance Lagrangian implicitly. 

The antisymmetric tensor formalism is used to describe the vector resonance since, as it is demonstrated 
in Ref. [49], only in this way one does not need to include extra terms in the pseudo-Goldstone Lagrangian 
to fulfill the QCD short distance constraints. For the vector resonances uj and we assume ideal mixing 
throughout 

= - Y^*/*, ^ \!\'^ ^ \[\^ ' ^^^^ 

and we do not include any 1/Nc suppressed operators in this respect. Nonetheless, we employ different 
bare masses for the p{770) and K* (892) in order to obtain a good fit to data. 

Before finishing this section, we introduce the last pieces of the chiral Lagrangian involving only 
pseudo-Goldstone bosons at 0{6) [40] 



C^ = A,—D^i;D^'i;-iA2 — {U+X-X^U)i; , 



^2 jp2 

—D^i;D^^i;-iA2 — 
■ip = —i Indet U , 

D^V = - 2(a^) , (12) 



with afj_ = {r^ — l^) /2. The inclusion of such operators does not improve our fits to data indeed. However, 
we take into account the monomial proportional to A2 to bring our prediction for the masses of r] and 
77' at their physical values |1J. This is necessary in the present work, since in the fit all of the pseudo- 
Goldstone masses have their physical values, while we need to use our prediction for the masses when 



^The term with Cm is purely proportional to rji while that with Cd requires at least one r;i, which mainly becomes an rj' 
because the r;i contribution to the rj is suppressed by the pseudoscalar mixing angle, as it is shown later. 
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we discuss the Nc dependence of our results, in particular the movement of resonance poles with Nc- 
Thus, it is necessary to match our prediction for the masses of the pseudoscalars as a function of Nc 
with their physical values for Nc = 3. Precisely the A2 term leads to an important contribution to 
the ry and r/' masses and in the triple expansion scheme its chiral order is lower than the one of chiral 
loops. Though this operator contributes to meson-meson scattering as well, it mainly contributes to the 
processes involving 77' and the inclusion of this term barely affects the global fit. Unlike the A2 operator 
the monomial proportional to Ai in Eq. (jl2p only contributes to the masses and the scattering amplitudes 
in an indirect way, i.e. through the normalization of the r/' field, and its influence in the global fit is tiny. 
Indeed, if we include this counterterm in our fits to data the resulting fitted value tends to vanish. As 
a result, we do not consider any further the Ai term neither in scattering nor for the rj — rj' mixing and 
masses. 

Finally, we want to point out that there is no double counting problem by having both the resonance 
Lagrangians in Eqs.dS])-® and the local pseudo-Goldstone operator A2, since the A2 term can be only 
generated by integrating out the excited pseudo scalar resonances, instead of the scalar and vector ones 
considered here. 

3 Structure of the scattering amplitudes 

Even in the leading order Lagrangian Eq.([T|), the fiavor eigenstates r/g and rji are not mass eigenstates 
and we use the angle 9 to describe the mixing of 773 and rji at this order 

r?8 = cefj + sgf]' , 

Vi = -sev + cel]', (13) 

with cg = cos 9 and sg = sin 9. In our notation rj and rj' are the fields that diagonalize the quadratic terms 
of the Lagrangian Eq. ([1]). The r] — rj' mixing at leading order is discussed in Appendix [Bl Eqs. (jB.SP - 
(lB.7p . The differences between r/,r/' and the physical states rj ,r/' are caused by higher order operators, 
including loops, and can be treated perturbatively within the triple expansion scheme. 

Next we calculate the contributions beyond the leading order to the scattering amplitudes in terms 
of the rj and rj' fields, while for the leading order terms to meson-meson scattering, stemming from Cy^ 
Eq. ([I]), one has to take care of the full r;— r/' mixing (see below). We also point out that for the calculation 
of the basic amplitudes it is more reasonable to use the and 77' fields than the 771 and 773 ones. This is 
because the insertion of the leading order mixing of 773 and 771 does not increase the order of a diagram, 
as illustrated in Fig. [1] for a one-loop contribution. The leading order mixing is proportional to m\ — m^, 
and is always accompanied by the inclusion of one extra 773 or 771 propagator that compensates the chiral 
power of the vertex. As a result loop diagrams with an arbitrary number of insertions of 773 — 7/1 mixing 
vertices have the same order. On the contrary, the mixing of 77 , 77' only receives contribution from higher 
orders, which guarantees that diagrams with insertions from the 77 — 77' mixing are indeed suppressed. 




Figure 1: The dot denotes the mixing of 773 and 771 at leading order, which is proportional to m?^^ — 
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The calculation of the scattering amplitudes comprises the contact vertex from £^ Eq.([T]), the one 
loop and resonance exchange graphs, as illustrated in Fig. [2j In addition, the mass and wave function 
renormalization terms, displayed in Fig. [3l should also be included. The latter ones only affect the tree 
level scattering amplitudes from Eq. ([T]) up to the order we attempt to calculate, as the contributions 
from the resonances and the loop diagrams are already beyond the leading order. The decay constants 
of the pseudo-Goldstone bosons, illustrated in Fig. HI also give rise to one higher order contributions 
by rewriting F in terms of the physical ones in the leading order scattering amplitudes. We express the 
scattering amplitudes employing one single physical decay constant F.,^ throughout. The relation between 
F and is given in Eq. (lC.2p . 





Figure 2: Relevant Feynman diagrams in the scattering amplitudes up to one-loop order. In digram (d) 
the coupling of the scalar resonances with the vacuum is indicated by a cross. 




Figure 3: Relevant Feynman diagrams for the pseudoscalar self-energy. 



S 



Figure 4: Relevant Feynman diagrams for the pseudoscalar decay constants. The wiggly line corresponds 
to the axial-vector external source. 



Another subtle contribution to the scattering amplitudes is related to the r] — r]' mixing, as commented 
above. We recall that rj ,rj' result from the diagonalization of r]s,r]i at leading order, Eq. (jl3p . After 
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including the higher order contributions from resonances and chiral loops, r] will mix again and the 
physical states r] , r( can be obtained by diagonalizing 77 , This extra diagonalization process contributes 
relevant pieces to the scattering amplitudes through the leading order results from Eq. ([TJ. Thus it is 
necessary to work out the rj— tJ' mixing at the one loop level by calculating the diagrams in Fig.[3l However, 
for the already next-to-leading order contributions and higher, namely, those amplitudes obtained from 
the exchange of resonances and involving loops, we do not need to distinguish between r] , rf and r/ , r/', as 
their differences only cause higher order effects that are beyond our current consideration. 
We parameterize the higher order rj — r]' mixing as 

^ -VV -(5^2 777/ . (14) 

where nijj and m^/ defined in Eqs. ()B.5P and ()B.6P stand for the leading order masses of rj and ry' respec- 
tively, while the different 5i, given in Eq. (IB.Sp . contain the higher order contributions. The physical 
eigenstates r/ and r/' are related with the rj and rj' fields by diagonalizing canonically the quadratic terms 
in Eq. (jl4p in the following way 



\ _ / cos Us — sm 
rj' J \ sin 9s cos 





(15) 



We calculate the leading order contributions, corresponding to the diagram (a) of Fig. [21 in terms of 
the physical r] and r/' fields, so that we can directly use the physical values for the masses. Up to the 
order considered in our calculation for meson- meson scattering, namely, up to 0{6^) for local tree level 
contributions and up to 0{6'^) for the one-loop graphs, we can invert Eq. (jl5p in a perturbative way. 



Then we have 



r?\_/l-% -f \ / cos^^ sm9s \fv 



2 



rj' J \ _^ I vy\— sin 6^ cos 9 s J \ v' 



(16) 



where 9s is determined through 



6„ 



,2 



tan6'5 = ^ —, (17) 

m„, — mp. 



with 



1 



rj' rj ' mf., rf V ) 



2m\ = rn\^ rn\, - ^ {iffi^ - m^^ y + 45^^ • (18) 

In the previous equations 9s and 5i are originated by the higher order contributions from the resonances 
and one loop graphs. Up to the precision we consider, it is safe to take only linear terms in the 6i so that 
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COS 6s = \/l — sin^ 9s — 1- From the relations between rj ,rj' and rj ,rj' , given in Eqs. (|15p and (jl6p . there 
are two ways to proceed for the calculation of the physical scattering amplitudes. One method, the one 
that we follow, consists of writing the physical amplitudes by expressing fj and fj' in terms of the physical 
fields T] and rj' in Eq.([T]), and then calculate the tree level amplitudes. In this way, all the masses from 
the kinematics are the physical ones, since they correspond to the rj and ij' fields. Another method is to 
determine the amplitudes with rj and rj' by using Eq. (ll5p as a linear superposition of those calculated in 
terms of the bar fields. However, one should bear in mind that for this case the masses of rj and rj' from 
the kinematics are those of the physical rj and r/' fields to which they are attached [45j. 

The explicit expressions for the mass and wave function renormalization, mixing parameters entering 
in Eq. p^ . pion decay constant, and scattering amplitudes are given in the Appendices [BtiDl 



4 Partial wave amplitude and its unitarization 

Once the perturbative amplitudes are calculated from U{3) xPT, as done in the previous section, we can 
proceed to perform the partial wave projections of the isospin amplitudes and construct the corresponding 
unitarized amplitudes as well. The amplitudes with well-defined isospin / from different processes are 
derived by assigning the following phase convention to the pseudoscalars in connection with the isospin 
basis states. This convention is consistent with the one taken in Eq. 

|r/) = |00), |r?') = |00), 

|7r+) = -|ll), |7r-) = |l -1), K°> = |10), 

\K-)^-\ll),\K^)^\l-'-)AK^)^\ll),\K-)^\l-l), (19) 

where I//3) is a state with isospin / and third component I3. 

For vrvr — t- vrvr scattering there are three isospin amplitudes, / = 0, 1, 2. They read 

T^{s,t,u) = 3A{s,t,u) + A{t,s,u) + A{u,t,s) , 
T^{s,t,u) = A{t, s,u) — A{u,t, s) , 

T'^{s,t,u) = A{t,s,u) + A{u,t,s) , (20) 

where A{s,t,u) stands for the process tt+tt > tt^tt^ and s,t,u are the standard Mandelstam variables. 

These equations, and other similar ones that follow, are obtained by invoking crossing symmetry |52j . 

For Kir — t- Kit and vrvr — t- KK, the different isospin amplitudes can also be expressed in terms of one 
single amplitude, again by using crossing symmetry. One then has 



Tl{s,t,u) 




+ {s,t,u) , 










2 


{u,t,s) 




T^{s,t,u) 




ri(t,s 


3 

u) +T2{u,s,t) 


T\s,t,u) 


= tI{v 


,s,t) - 


Tl{t,s,u). 


For vrvr — ?■ 7777 and vrr/ — )• vrr/, we have 









(21) 



r°(s,t,M) = -V3C{s,t,u), 

T^{s,t,u) = C{t,s,u), (22) 
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with C{s,t,u) = T^o^o^j^fj{s,t,u). For vrvr — )■ r]r]',r]'r]' and vrr/(ry') — )■ vr?]', analogous formulas emerge and 
the only difference is that C{s,t,u) is redefined accordingly to the scattering process. 
The reactions Kn — )• Krj and Krj' are pure 1 = 1/2 and are given by 

{s, t, u) = -V3Tk+^o^k+,j{s, t, u) (23) 

and similarly when the 77' is in the final state. In terms of them one also has the scattering amplitudes 
for the t-crossed processes KK — )• ttt] and vrr/', which is pure 1=1. It reads 

T\s, t, u) = -V2TK+^o^K+^{t, s, u) , (24) 

and analogously for KK — )• nrj'. 

There are two isospin amplitudes in KK — t- KK, that can be expressed as 

r°(s,t,M) = 2D{s,t,u) + D{t,s,u) , 

T^{s,t,u) = D{t,s,u), (25) 

with D{s, t, u) = Tj^+K-^KORoi^j t, u). Notice that here we have derived both isospin amplitudes in terms 
of only one, while previous works used both r^+^_^^o/fo(s, t, w) and Tj^+j^-^j^+f^-{s,t,u) [53 1 [51 1 [55] . 
For KK —7- rjr] and Ki] — t- Kj], the amplitudes involved read 

T°(s,t,u) = -V2E{s,t,u) , 

T^{s,t,u) = E{t,s,u). (26) 

with E{s,t,u) = Tj^oj^o^,^^{s,t,u). For KK — )• rjrj',r]'r]' and Kr]{rj') — Krj and Krj', the formulas are 
analogous and the only difference is that one needs to redefine E{s, t, u) for the corresponding process. 

All the processes r]^'^ri^'^ — )• 7]^'^r]^'\ with ry^'-* either an rj or an r/', are / = and no crossing relation 
can be invoked to reduce the number of amplitudes needed, one for each process. 

Next we can perform the partial wave projection for the isospin amplitudes and the convention we use 

is 

Tjis) = l\xPj{x)T'{s,t{x),u{x)) , (27) 

where x = cos ip, with if the scattering angle between the incoming and outgoing particles in the CM, 
and Pj{x) stands for the Legendre polynomial. In the previous equation (-v/2)^ is a symmetry factor to 
account for identical particle states, such as vrvr (with isospin 0, 1 or 2), r/r/ and r/'r/'. It corresponds to 
the so-called unitary normalization of Ref. [9j. E.g. = 1 for vrvr — )• KK and N = 2 for vrvr — )■ r/r/, and 
so on. 

The unitarization method we follow was developed in Ref. [IQJ and is based on the N/D method [llj . 
The essential of this approach is to separate the crossed-channel cuts (LHC) and the right-hand cut 
(RHC), due to unitarity, in two different functions, N{s) and g{s). The former does not contain the two- 
particle unitarity cut but it has the LHC, while the later contains only the two-particle unitarity cut and 
not the LHC. A more detailed account of this unitarization method can be found in Refs. [10 1l56[[571 l58j. 
In the latter reference an explicit integral equation for the function N{s) in nucleon-nucleon scattering 
is deduced. 

Let us consider first the elastic case, its generalization to the coupled channel case is straightforward 
and given below. Because of unitarity above the two-particle threshold sth a two-meson partial wave 
amplitude Tj{s), with well-defined isospin / and angular momentum J, fulfills 

Imri(5) = T'j{s) p{s) ri(s)* ^ Imri(.)-i = -p(s) . (28) 



10 



In the previous equation 



, . \/\s - {ma + nibYWs - {nia - mfe)^] 



167rs 

(29) 



with nia and m?, the masses of the two particles in the state and q the CM three-momentum. Eq. (|28p 
imphes that the imaginary part of the inverse of a partial wave is known. This can be used to write down 
a dispersion relation of the inverse of Tj{s) taking as integration contour a circle of infinity radius that 
engulfs the right hand cut jlO) . It then results 

Tj{s)-' = N^is)-^ + g{s) , (30) 

so that the function of Nj{s), by construction, does not contain the RHC and g{s) results from the known 
discontinuity along this cut, Eq. (|28p . In this way g{s) is given by the following dispersion relation, 

g{s) = g{so) - ^ r P[f] -ds' . (31) 

^ ^^th [s' - s){s' - So) 

with 

Imgis) = -pis) , (32) 

for s > Sth = {ma + m{,)^. The dispersive integral in Eq. ([5T]) can be represented by the typical two-point 
one loop function |10] 

x± = ^ + "^b ^ J_^-4s{ml - iO+) + {s + ml- mg)^ , 
167r^5(s) = asL(/i) -Flog ^ -x+ log x_ log — . (33) 

X-\- X- 

For the case of equal mass scattering, g{s) reduces to the simple form 

^(^) - 1 
/x^ '''''''' a{s) + I 

with 



levr^ g{s) = asUp) + log — - ^{s) log J} J, , ^ , (34) 



a{s) = ^l-^. (35) 

In order to determine the interacting kernel Nj{s) we proceed as explained in Refs. |56l I57j . so that 
we match the unitarized amplitude Eq. ()30p to the perturbative calculation, which is performed up to the 
one-loop levelJl In this way, Eq. ([30|) should be expanded up to one power of g{s). It then results for 

nK^), 

iVj(s) = rj(s)(2)+Rc«+Loop + rj(s)(2) g{s) Tj{s)^^^ . (36) 



^Eq. (I30|l is valid below the first threshold of multi-particle states. Note that such states are further suppressed in the S 
counting beyond our present one-loop calculation for the interaction kernel (they at least imply a two-loop calculation) and 
we do not consider them any longer. 
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Here rj(s)(2)+R^'^+Loop stands for the partial wave amplitude from the perturbative calculation, with the 
superscripts (2), Res and Loop denoting the tree level amplitude from Eq. ([1]), resonance exchanges and 
loop contributions, respectively, depicted in Fig.[2j The wave-function renormalization and r] — r]' mixing 
contributions are included in Res and Loop. 

It can be easily checked that Nj{s) from Eq. ()36p does not contain the RHC. This is due to the fact 
that the perturbative partial wave amplitudes satisfy unitarity perturbatively 

imTi(s)(2)+R-+L°°p = ri(s)(2) ri(s)(2) 

= -TM'^lmgis)Tjis)('\ (37) 

for s > Sth- Regarding the LHC contribution, we would like to emphasize that it is perturbatively 
treated and collected in the Nj{s) function defined in Eq. (l36p . unlike the RHC that we take into account 
non-perturbatively. 

If the amplitude Tj{s) has a zero in the complex plane, e.g. an Adler zero that appears in the 
subthreshold region of some partial S- waves as a consequence of chiral symmetry [59) . this is accounted 
for in Eq. (I30p by the corresponding zero of Nj{s). The zeroes on the real axis correspond to poles of 
the inverse of the amplitude, the so-called Castillejo-Dalitz-Dyson poles |60j . As discussed in detail in 
Ref. |T0] when the LHC contributions are neglected Tj{s) = Tj{s) / (f"^ {s) can be written as 



ms) - -^^^ r ds'^^q^M^ . E a..- + i: ^ . (38) 

■J Sth ^ m=o i * 

In the previous equation J + 1 subtractions at sq have been taken, being the Om, m = 1, . . . , J the 
corresponding subtraction constants. The last sum accounts for the presence of poles in the inverse of 
the partial wave, with Si and Ri their locations and residues, respectively. For the S-waves one has the 
location of non-trivial poles due to the Adler zeroes. For a given partial wave its position fixes one of 
the Si pole positions in the right-hand-side of Eq. ([55]) . In the presence of LHC, Dj{s) = Nj{s)~^ ~^ di^) 
by comparing the previous equation with Eq. (j30p . the zero in Nj{s) is moved according to the chiral 
expansion, Eq. (j36p . from its leading position given by Tj{sY'^\ In the case of the lAM one has to modify 
the standard formula [5l [U [7] to account properly for the Adler zero region [61] , otherwise one has a 
double zero instead of a simple Adler zero and a spurious pole in the partial wave amplitude in the same 
subthreshold region [62] . 

The previous formalism can be easily generalized to scattering processes with multiple coupled channels 
by employing a matrix notation. In this way, Tj{s), Nj{s) and g{s) — gj(s) are now matrices and Eq. ([30]) 
still holds. Since phase space is diagonal then the matrix gj{s) is also diagonal, with its matrix elements 
given by Eq. ()33p . evaluated with the appropriate masses for the corresponding channel. The matrix 
Tj(s) is symmetric, as required by time reversal invariance [3], which then implies from Eq. (j30p that 
Nj(s) is also symmetric. 

Explicitly, for the IJ = 00 case there are five channels and the corresponding matrices read 



( ^-K-K^-K-K ^ -K-K^KK ^-KTr^riT] N-KTT^rjr)' ^TT-K^ri'ri' ^ 

^TTTT^KK ^KK^KK ^KK^rirj ^KK-^qr)' ^KK^r)'r)' 

^TTTT-^rirj KK^rjT] -^rjri—^rjrj ^rj-q—^rirj' -^rirj—^rj' ri' 

-^TT-K^'qrj' ^KK—^rjrj' ^•qri^Tjri' ^rjr]' ^'qrj' ^'qrj'-^ri'ri' 



(39) 
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For the / J = 1 channel, the matrices are 



/ 


9iTn 














\ 







9kk 





















drj-q 





















9r]T]' 







V 
















/ 



N. 



N. 



N. 



^nri^KK ^KK^KK ^KK-^irri' 



N. 



7r?7'— >7r?7' 














9kk 











9nr]' 



For the / J = i channel, it results 



Nk-k^Kt^ Nk-k-^Kt] NKTT^Krf 
Nk-k^Kt) Nkti^Kt] Nx-q^Kri' 
Nktt^Kt]' Nkti^Kti' Nkti'^Kt]' 



dKir 











9Kn 











9Kri' 



Similar results hold for the IJ 



5 1 channel. 
For the IJ = 11 quantum numbers we have 



TTTT—^KK 



9l{s) 



?7r7r 

gKK 



(40) 



(41) 



(42) 



(43) 



(44) 



(45) 



(46) 



Let us consider now the purely elastic channels where Nj{s) and 9j{s) are just functions, not matrices. 
)r / J = I 0, one has 

nI{s) =NKn^K^, (47) 

3 

5o(s) =9K7r- (48) 

For IJ = 2 0, 



-^0 (*) ~ -^TTTT-^-TTTr ) 
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9o(s) =9nn- (49) 

Finally, the appropriate functions for JJ = 1 are 



9\{s)=9kk- (50) 

After having the unitarized scattering amplitude from Eq. (|30p . the S-matrix for the IJ channel, 5'j(s), 
can be defined straightforwardly in matrix notation 



55 = 1 + H^p'j{s) ■ Tj{s) ■ ^pijis) . (51) 

with yOj(s) = —\mg^j{s). From the matrix elements of the S-matrix we can read out the phase shifts 5kk 
and 5ki, with k ^ I, since 

Skk = I'S'fcfcle^*^'''' , 

Ski = ISkile'^" ■ (52) 



5 Discussion and results 

In this section we first discuss the fit to experimental data in order to fix the free parameters in our 
approach. Later we discuss the associated spectroscopy and its properties under variation of Nq- 



5.1 Fit Quality 

We perform the fit for the IJ = 00 channel up to y/s = 1300 MeV. The inclusion of the rj' meson is not 
enough to guarantee that our current calculation can be applied higher in the energy region. There are 
good phenomenological reasons to expect that the 4tt state plays an important at this energy level and 
its infiuence can not be simply neglected [12]. The observables fitted are shown in Fig. [5] and correspond 
to the elastic tttt — )• vrvr phase shifts, J^J'^^tttd the elasticity parameter, IS*^^.^^^!, and the phase and one 
half of the modulus of the S-matrix element for the inelastic process vrvr — )• KK above the KK threshold, 
^^^T^^KK ^^'^ ll'^^^-s>K_ff I' ^sspectively. For references to the experimental data see Fig. [5l 
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Figure 5: (Color online. 

phase shifts of irn — >• vrvr the modulus of the S-matrix element for vrvr — )• tttt (|S^^_5.^^ 

half modulus of the S-matrix element for tttt — )• KK (|5^^_^^j^|/2) and the phase shifts of vrvr - 
(5™ rt-E?)- ^JJ^-s-TTTT data correspond to Ref. [63] (triangle in green), [6l] (square in blue) and the average 



00 case. From top to bottom and left to right: the 

the 
KK 



00 



is from Ref. 



data from Refs. [Ml [Ml [67] (circle in black), as employed in Ref. [TO]. \S[ 
|5*"^ „i>|/2 is from Refs. [68] (square in blue) and [69j (circle in black). The phase shifts r^f> are 



from Refs. [M] (square in blue), [70] (circle in black). The solid (red) line corresponds to the best fit, 
Eq. (jSSp . while the error bands are presented by the shadowed area. Finally, the more constrained fit of 
Eq. (|59p is given by the dashed line. This notation also applies in Figs. [6] and [3 



For the IJ = ^ channel, since there is no significant inelasticity above the Krj' threshold, we fit the 
data up to y/s = 1600 MeV. Although for higher energies one already has the influence of the /Cg(1950) 
resonance [13] which we have not included. The observables in the fit are the Kir — )• i^vr phase shifts 
io 

^K-K^Kn- '^^^ reference to experimental data is given in Fig. [6l 

For the I, J =10 channel, there are no available data for scattering up to now and it is typical to use 
a vr?7 event distribution from other production processes to obtain information on the ao(980) resonance. 
As the production mechanism is not under good theoretical control in our work, we decide to fit the 
data around the resonance region up to 1050 MeV, as explained below. The experiment data for the irr] 
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distribution are illustrated in Fig. [6l 




Figure 6: (Color online.) From top to bottom and left to right: phase shifts of Ktt — )■ Ktt with IJ = 2 ^ 

1q 

{^k-k^k-k)i event distribution with IJ = 10, phase shifts of vrvr — )• vrvr with IJ = 20 and 

— —0 

phase shifts of Ktt — )■ Ki: with / J = | {^K-K^Kn)- ^k-k^k-k corresponds to the average data from 
Refs. [Til [72l [73] (square in blue), as used in Ref. [10], and Ref.[74] (circle in black). Data for the vrr/ 
event distribution are from Ref. [75] and the dotted line corresponds to the background [10]. 5"^^^^^ is 

from Refs. [76] (square in blue) and [77] (circle in black). The experimental data for Sf^^^j^^ are taken 
from Refs. [75] (square in green), |79j (circle in blue) and [72] (triangle in black). For the meaning of 
lines see Fig. [5l 



For the exotic channels with I J = | and / J = 2 0, we fit the phase shifts of Ktt — )• Ktt and vrvr — t- vrvr, 
respectively. The references for the experimental data are given in Fig. [6l 

The observables we fit in IJ = 1 1 and U = \ l are the phase shifts of vrvr — )• vrvr and Ktt — )■ Ktt, 

5}^^^^^ and ^K-K^K-K^ order. They are fitted up to 1200 MeV, as we only include the lowest multiplet of 
vector resonances in the Lagrangian Eq.([6|). The references to experimental data can be found in Fig. [71 
Before moving to the details of the phenomenological discussion, let us first comment on the free 
parameters in our theory. Apart from the couplings in Eqs.([5]), ([6|) and (jl2p . there are the bare mass 
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parameters of the resonances, the U^i^) anomaly mass Mq in Eq.([T|) and the subtraction constants asL 
defined in Eq.(|33p for different channels. One way to reduce the number of the subtraction constants 
is to impose the isospin symmetry for them [80J. In this way, 05^ of ^tttt in I J = 00 and in I, J = 20 
are the same. Similarly, Qxk = 00 and IJ = 10 are also equal. This also applies to qk-k in 

IJ = ^ and in IJ = | 0. In principle, the subtraction constants in the remaining channels are free. 
Nevertheless, if one further assumes ?7(3) symmetry within each orbital angular momentum the values of 
different subtraction constants should be also equal. For this result one should extend the proof given in 
Appendix A of Ref. [80] in the SU{'i) limit to the C/(3) case, that holds for equal quark masses at leading 
order in large Nc- Notice that in this limit the quarks behave equivalently so that the equality of the 
subtraction constants in SU{2>) also implies that they are equal for the U{2>) case. We indeed exploit this 
feature as much as we can in the numerical discussion. We summarize the subtraction constants used in 
the fits discussed below: 
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(53) 



and all of the subtraction constants in the vector channels are set equal to a^^. 
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Figure 7: (Color online.) Plots for the vector channels from the fits, Eqs. (I55p and (I59p . The left panel 
11 channel and the right one is for IJ = ^1- The data of phase shifts for vrvr — ?■ vrvr with 
are from Refs. [81] (square in blue) and [82] (circle in black). The data of phase shifts 

are from Refs. (circle in black) and ^72j (square in blue). For 



is for the I J 
I J = 11, J" 



TTTT — f-irn 



for Kir — )• Kir with I J 



'^Kn- 



the meaning of lines see Fig. [5l 



In order to fit the ttt] mass distribution we need to introduce two additional parameters to parameterize 
the production mechanism, J\f and c, that enter in the expression 

= Qnv^l TKK^nvis) + cT^r^^^r,{s) f ■ (54) 
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Here, Qt^y) is the three momentum of ttt] system in CM and E-j^r) = \/s is the energy in the same frame. The 
parameter M accounts for the fact that the event distribution is not normahzed. The hnear combination 
of amphtudes in the previous equation originates because the ao(980) shows up differently in them. See 
also Ref. [56] for a more detailed explanation of the derivation of Eq. (I54p in connection with the invariant 
mass distribution of KN around the A(1405) resonance. 

In the end, we have 16 free parameters and the fitted results are 



with xVd.o.f = 714/(348 



Cd = 


(15.6i^:^)MeV, 


Cm — 


(31.5ti9|)MeV, 


Cd = 


(8.7l?:^) MeV , 


Cm - 


= (15.8l|3)MeV, 


Mss 


= (1370t^f )MeV, 


Ms, 


= (1063115) MeV, 


Mp-- 


= (801.01^:^) MeV, 


Mk* 


= (909.0l^j)MeV 


Gv- 


= (61.9l}|)MeV, 


10,7rr) _ r, r) + 3.1 

"5L ~ ^•'-'-3.4 ' 




= (-1-151^9'), 


<^SL 


= (-o.96i°:}°) , 


M = 


(0.6l|]J) MeV-2 , 


c=( 




Mo-- 


= (954;^f )MeV, 


A2 = 


= (-o.6i°i) , 


-16) 


~ 2.15. 







The corresponding figures from the fit are displayed in Figs. O [6] and [71 The width of the bands shown 
and the errors given in the fitted parameters in Eq. (j55p represent our statistical uncertainties at the 
level of two standard deviations [70] 

= AxV(2xg)'/' , (56) 
being Xo the minimum of the obtained corresponding to the best fit, the number of standard 
deviations and Ax^ = ~ Xo- 

As one can see in Figs. El [6] and [71 most of the observables are well reproduced through the fit. The 
phase shifts of vrvr — )■ vrvr with IJ = 00 in the low energy region from the recent measurement of Kg^i^ decay 
|64j are included in our fit and they can be perfectly reproduced. For the Trry invariant mass distribution, 
we explicitly subtract the background as analyzed in Ref. [75], which are mainly caused by the tail of 
higher resonances. Compared with observables in the scalar channels, the errors in the vector channels 
are rather small. Among all of the curves, the least satisfactory results from our fit correspond to the 
channels with exotic quantum numbers with IJ = 20 and IJ = |0, which is also an important source for 
the x^ resulting from our fit. For example, with 42 data points, they contribute 281 to the total x^j i-e. 
12% of degrees of freedom carry 39% of the x^- Nevertheless, one also has to remark that the different 
experimental data in these channels are not always compatible at the level of one a, which also makes the 
X^ to increase. On the other hand, since no resonances appear in those exotic channels, it may indicate 
that unitarity is less important as compared with the other channels. In contrast, the LHC contribution 
could play a more important role in this case |101 1^ . We recall that precisely our unitarization scheme 
only incorporates the LHC effects perturbatively. Thus, an improvement on the treatment of the LHC 
should provide a better fit for these two channels. 

With the fitted parameters in Eq. (I55p we can also calculate the values for the masses of the r] and ry' 
that are obtained by diagonalizing Eq. (|14p . We then obtain from Eq. (|18p the values 

= (526^^5) MeV , m^, = (951^^^) MeV . (57) 

The leading order mixing angle of r/ — r/' that results from Eq. dH) is 

e = -l%.2°tlT° . (58) 
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Concerning the resonance parameter Cd in Eq. (j55p . the value from our fit is a bit smaller than the 
previous determinations, such as q = 19.1^2'^ MeV given in Ref. [10], q = 23.8 MeV reported in 
Ref. [l5] and Cd = (26 ± 7) MeV in Ref. [83]. Nevertheless, they are compatible within the error bands. 
For the coupling Cm, though different approaches predict a broad range for its central value, they always 
accompany a large error, for example Cm = 31.5^22 5 MeV in the present case, Cm = (15 ± 30) MeV in 
Ref. [To] and Cm = (80 ±20) MeV given in Ref. [83]. Comparing the couplings related to the singlet scalar 
with Ref. do], our determination of Cm is compatible with Cm = 10.6l3;5 MeV in Ref. [TO], while for Cd 
the current fit leads to a smaller value, compared with Cd = 20.9lj^'Q MeV given in the same reference. 
For the bare masses of resonances, our present results agree well with the previous determinations: 
Msg = (1390 ± 20) MeV, Ms^ = (1020t^[j) MeV in Ref. [lOj and Ms^ = 1400 MeV in Ref. [l5]. Our 
values in Eq. (f55]l are compatible with the estimates Cd — Cdl\/'i> and Cm — Cm/V^ based on large Nc 
|48j . however the bare masses for the singlet and octet scalar resonances are different, at around a 30%, 
as obtained in previous studies [101 112j . The current determination for the subtraction constant asL is 
consistent with the previous result in Ref.[10| as well. For the vector resonance coupling Gy-, our result 
is in good agreement with Gy = (63.9 it 0.6) MeV determined in Ref. [83j. The C^a(1) anomaly mass Mq 
carries a large error and agrees with the conclusion in Ref. [33]. The 1/Nc suppressed parameter A2 is 
poorly known in the literature and in Ref. [^ A2 = 0.3 is estimated, which is somewhat incompatible 
with ours. Nevertheless in Ref. |47j the value of A2 is not determined by any physical observables, but 
only by naive dimensional analysis. Moreover in their case the counterterm A2 is always accompanied by 
the factor of m^, indicating the insensitivity of this parameter there. 

We also show another more constrained fit by employing relations and values in the literature, already 
commented above, for the parameters shown in the fit of Eq. ()55p . In particular we impose from the 
beginning the large Nc constraints [31] q = Cd/^/S and Cm = Cm/V^- We also take M^g = 1390 MeV 
and = 1020 MeV, according to Ref. [TU] . For Gy we take 60 MeV from the averages of values taken 
from Refs. [33 [23l [83]. For Mq we take the value 850 MeV from Ref. [3i]. The subtraction constant 
CL^gi^^ , given with large errors in Eq. (j55p . is now fixed at +2. The resulting fit is shown by the dashed 
lines in Figs. [5]|71 having now only 9 free parameters. We see that the reproduction of data is of similar 
quality as the one achieved by the fit in Eq. (|55]l and the x^/d.o.f = 843/(348 — 9) = 2.5 is not much 
larger. The resulting values for the free parameters are now 

Cd = 17 A MeV Cm = 28.1 MeV 
Mp = 800.4 MeV Mk* = 910.0 MeV 

af^ = -1.14 all = -0.89 (59) 

A2 = -0.22 Af = 0.55 MeV"^ 

c = 0.84 

Comparing with the values in Eq. (|55p one observes rather similar values, compatible within the errors 
given in Eq. (j55p . The biggest change occurs for the value of A2, although the new value lies well within 
the error band given in Eq. (1551) . The resulting singlet couplings are Cd = Cd/VS = 10.05 MeV and 
Cm = Cm/V3 = 16.20 MeV. 

5.2 Resonances generated from unitarized amplitudes 

We summarize the masses, widths and residues of the various resonances in Table [H Resonances are 
characterized by their pole positions in the partial wave amplitudes in unphysical Riemann sheets. Around 
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Table 1: Pole positions for the different resonances in = (M,-iE). The mass (M) and the half width (r/2) are 
given in units of MeV. The modulus of the square root of a residue is given in units of GeV, which corresponds to 
the coupling of the resonance with the first channel (specified inside the parentheses) . The last two columns are the 
ratios of the coupling strengths of the same resonance to the remaining channels with respect to the first one. The 
corresponding Riemann sheets where the resonance poles are located are explained in detail in the text. Note that 
here the residues for tttt, r]r] and rj'r]' are given in the unitary normalization, due to the extra factors of dividing 
Eq. (P?]). Thus, one should multiply by y/2 these couplings if one wishes to restore standard physical normalization 
to 1 for these states. The error bands of the resonance parameters appearing in the table only correspond to the 
statistical error from the fit in Eq. ([55|) . 



a resonance pole sr, corresponding to a resonance R, the partial wave amplitude Tj{s)i-^j tends to 

s - sr 

By calculating the residue of the resonance pole we then obtain the product of the couplings to the 
corresponding decay modes. At the practical level we calculate the residues by applying the Cauchy 
integral formula 

gR^iQR^j = -TT-- f Tj{s)i^jds . (61) 

For every pole one has further to indicate in which unphysical Riemann sheet it lies. Each function g{s) 
has 2 sheets and their relation is given by [9] 

gii{s) = gi{s) + 2ip{s) , (62) 

with gii{s) the function analytically extrapolated to its second Riemann sheet and gi{s) the function in its 
first Riemann sheet, given in Eq. (j33p . Different Riemann sheets are easily accessed by deciding on which 
sheet every gi{s) function, associated to channel i, is calculated. In this way, for an I, J channel there are 
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2 possible sheets, with N the number of coupled states with the same IJ quantum numbers. Along the 
real s-axis above threshold changing sheet implies to reverse the sign of the imaginary part of the gi{s) 
function. In the following we conventionally label the physical or first Riemann sheet as (+, +,+,.. .). 
The second Riemann sheet can be reached by changing the sign of the first momentum, which is labeled 
as (—,+,+,+,...). The third, fourth and fifth sheets, also considered in this work, correspond to 
(— , — , +,+,.. .), (+, — , +,+,.. .) and (— , — , —,+,.. .), in order. More sheets can be obtained by taking 
more combinations of plus and minus signs between the brackets. At a given energy value y^, there is 
one unphysical sheet to which one can directly access from the physical sheet by crossing from s + iO^ to 
s — iO'^ the branch cut between the two thresholds T„ and T^+i, with T„ < ^/s < Tn+i- In our current 
notation, this specific unphysical sheet corresponds to changing the signs of all the three-momenta below 
the considered energy point ^/s, i.e 

(-^-^, +,+,...) (63) 

n 

In fact this is also the most relevant Riemann sheet where one should search the pole for a given resonance 
affecting that energy region, although other shadow poles may also appear in other unphysical sheets |84j . 
Thus in the following, we mainly present our findings for a given resonance on the complex plane in this 
specific sheet and simply comment the results appearing in other unphysical sheets. 

Through the unitarization procedure described previously, we can simultaneously get the relevant 
resonances in the considered energy region, such as cr, /o(980) and /o(1370) in the IJ = 00 channel, k 
and ^(5(1430) in /J = i 0, ao(980) and ao(1450) in / J = 1 0, 0(1020) in IJ = 01, p(770) in I J = 11 and 
K*{892) with IJ=^1. The pole positions and residues of these resonances are collected in Tabled! We 
comment on them channel by channel next. 

I J = 00 

Three kinds of resonance poles in the complex plane have been found, which we identify with the a, /o(980) 
and /o(1370) resonances. According to our previous discussion, the most relevant Riemann sheets for 

/o(980) and /o(1370) correspond to those with (— , +, +, +, +), (— , +, +, +, +) and (— ,—,—,+, +), 
respectively. Their pole positions are compatible within errors with those given in the PDG [1]. The pole 
for the a resonance in Tabled] is in close agreement with the pole position (484 it 14 — 255 it 10 i) MeV 
of Ref. [85], and slightly narrower, considering errors, than that of Ref. [13] (441]']g^ — 272^2^ i) MeV. 

From the residues given in Table [H one can see that the biggest coupling of the a resonance is to the 
vrvr channel. Nonetheless, it also has a large coupling to the KK channel. From the value in Tabled] one 
has 

compatible with 0.37ib0.06 obtained in Ref. [86]. In this reference this value is interpreted as an indication 
that the a resonance has a strong glueball component. 

One observes from Table [1] that the a couples very weakly to r]r]. This is typical of unitarized xPT 
studies |87j . However, from QCD spectral sum rules [881 189] the presence around 1 GeV of a broad glue- 
ball, (Tb, together with a, 82 = {uu + dd)/y/2 quarkonium, is required in order to fulfill the corresponding 
scalar sum rules. If the a resonance were purely gluonium then it would be an SU{?>) singlet and its 
coupling to pseudoscalar pairs of vr, K and r/g would be universal. In this case the coupling of the a to 
•qr] would be the same as that to vr^vr^. On the other hand, in Ref. [89j a maximal mixing between gb 
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and S2 is proposed, with a mixing angle 9s — (40 — 45)°. In this case the size of the couphng of the a to 
rjrj depends strongly on the value taken for the coupling of the a to a kaon pair. E.g. if one takes as in 
Rcf. [89] that \gcr^TT+TT- /9<7^k+k- \ ^ li foi^ definiteness we employ 1.2, then 



, ~ 0.8 . 

However, if the latest value 0.37 [86] for \ga^K+ K- / 9a^iT+n- 



is considered then 



0.15 . 



This result is remarkably close to our prediction from Table [T] 
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In the previous equation, as well as in the following, a tilde over a coupling of the a to a state made of 
identical mesons means that it is considered with physical normalization and not with the unitary one, 
used in Table [H see Eq. 07D . 

Up to the authors' knowledge this is the first time in the literature that the couplings of the a to states 
involving the pseudoscalar singlet rji are calculated. In Ref. [5H] one has the upper limit 
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< 0.23 . 



(68) 



From this value one can obtain the corresponding upper limit for the coupling of the a to r/'r/' dividing 
by ianO the one to r/r/'H Taking for 9 our value given in Eq. (j58p we then obtain the upper bound 
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From the values obtained in Table [J we have 
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which are perfectly compatible with the bounds given in Eqs. ()68p and ()69p calculated from ef. |88) . 

We then obtain a remarkable compatibility between our results for the couplings of the a in table [1] 
and those obtained from QCD spectral sum rules [88^ I89j. once the latest value for the ratio of couplings 
\9a^iT+iT- /9a^K+K-\ is taken from Ref.[86j. Nevertheless, in our approach there is no a corresponding 
bare state at the Lagrangian from which the a pole could stem, it is mainly generated by the strong 
interactions between pseudo-Goldstone bosons. 

The /o(980) resonance has its strongest coupling to KK, though it also couples almost as equally 
strong to r]r]. Notice that the ratios of the couplings to KK and vrvr between the a and /o(980) are 
nearly inverse each other. The importance of the rji] channel for understanding properly the /o(980) was 
stressed in Ref. [10], because once this channel is considered the inclusion of a bare single state around 
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1 GeV is required [121 HZ]- The /o(980) pole in fact moves continuously to the bare state around 
1 GeV if the (7o(s) matrix is removed, e.g. by multiplying it with a factor A G [0, 1] and taking A — t- 0. 
But let us remark that the strong KK interactions give rise to a bound state close to the actual mass of 
this resonance P [Ml EH ESI [93] , around the KK threshold. 

The /o(1370) resonance owns its origin to the octet scalar erg with the bare mass of 1370 MeV, as 
already noticed in Ref. [12]. These bare poles gain their widths through the unitarization procedure. The 
/o(1370) couples most strongly with the tjt]' channel. However, a deeper study of the /o(1370) resonance 
requires to include explicitly the 47r channel [12]. 

These conclusions are further supported by the Nq trajectories of those poles, which we will discuss 
later. 

/J=iO 

For these quantum numbers we obtain two poles corresponding to the k and -frQ(1430) resonances. 
The most relevant Riemann sheets correspond to (—,+,+) and (—,—,+) respectively, with poles at 
^ = (665+^ - 268lf i) MeV and ^/s = (1428^^^ - 87t^| z) MeV, in order. Their pole positions are 
compatible with those in the PDG [T] . For the k resonance the pole position (658 ± 13 — 279 it 12 i) MeV 
is found in Ref. [15] from a Roy-Steiner representation of Kir scattering, also in good agreement with our 
result. In the fourth Riemann sheet, i.e with (+, — , +), a shadow pole at ^/s = (717^5 — 280^24 i) MeV 
is also found for the k resonance. This resonance couples most strongly with the Ktt channel, although 
it has also a large coupling to the other channels. In turn, the ii'Q(1430) has its largest coupling to the 
Kr]' channel, which has its threshold quite close to the resonance mass. 

We also find a good agreement with the previous study of Ref. [45] that shares for this channel many 
facts in common with ours. The k resonance, similar to the case of a, is generated mainly from the 
pseudo-Goldstone interactions. The -fCQ(1430) pole originates from the bare octet of scalar resonances in 
the Lagrangian Eq. ([5|). Again, these conclusions are further supported below when discussing the Nq 
dependence of the pole positions. 

/J = 10 

This is the most problematic channel in our analysis. By imposing the U (3) symmetry on the subtraction 
constants, i.e. with the constraint of a]?/^^ = Og'l^^ = cls'l^^ ' good fit for this channel could be 
obtained simultaneously with the other data. Within a reasonable range for the ttt] subtraction constant 
the best fit prefers positive values, in contrast with the negative ones for the other channels, though the 
error for this parameter is large. However, we notice that, similarly to previous studies [9l [lOj, if only 
the tree level amplitudes are taken for the interacting kernel one can fit well the same data and get a 
reasonable pole for the ao(980) resonance on the second sheet. Also a similar negative value for the irr] 
subtraction constant results then as in the other channels, with the constraint a]?l^^ = = Og'i^^ 

|10] . We have checked too that the influence of the ttt]' channel to the mass distribution of vrr/ Eq. ()54p is 
actually tiny. In addition, no two-body scattering data are available for this channel which also reduces 
its statistical weight in the fits performed. 

We show in Table [1] the ao(980) and ao(1450) pole positions and residues. For ao(980), the most 
relevant sheet should correspond to the one with (— , +, +) or (— , — , +). However with our current best 
fit, no reasonable poles are found in those Riemann sheets. Nevertheless, a pole located in the fourth 
sheet, i.e. with (+, -, +), at ^ = (1012+^ _ iQ+fO i) MeV appears. This rather indirect location of the 
pole associated with the ao(980) resonance indicates that this resonance is to a large extent a dynamical 
cusp effect within our approach. It couples mostly with the KK channel as shown in Table [TJ 
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The most relevant Riemann sheet for the ao(1450) corresponds to the one with (—,—,—), since its 
mass is larger than the Trry' threshold, where one has the pole ^/s = (ISGS+gg - Jl^H i). Other channels 
are lacking at the typical energies for the ao(1450), as listed in the PDG [1] for its decay widths. This 
lack of decay channels (e.g. ojutt and ao(980)-7r7r) explains why our pole position for this resonance gives 
a smaller width than in [1], though the mass is compatible within errors. In our approach this resonance 
dues its origin to the bare octet of scalar resonances in the Lagrangian of Eq. ([5]) . 

IJ = 11, \ l and 01 

For the vector channels we find poles corresponding to the p(770) and iC*(892) in their respective 2nd 
Riemann sheets at = (762^^ - 72^2 i) MeV and ^/s = (891+^ - 25+^ i) MeV, respectively. Data are 
reproduced in a straightforward way, though we have to distinguish between the bare p and K* masses 
with fitted values given in Eq. (I55p . The quality of the fit is the same regardless we free their subtraction 
constants or fix them to the scalar channels. This indicates that unitarity effects are not so important, 
though they provide the right widths to the bare poles [10]. The masses and widths from the resonance 
poles agree well with the PDG values [1]. 

In addition, we have another pole corresponding to the </>(1020) resonance in the vector isoscalar 
channel. We obtain a width for this resonance to two kaons of around 4 MeV, which is close to the 
experimental partial decay width of the (/)(1020) to this decay channel of 3.5 MeV pLj. Nonetheless, 
another important decay channel for the (/)(1020) is the Trvrvr state which is not considered in our approach. 



5.3 Nc trajectories of the resonance poles 

As discussed above in the Introduction, one of the important improvements of our current work is to 
take into account the Nq dependence of the pseudo-Goldstone boson masses when discussing the Nc 
trajectories for the resonances. This may potentially cause some significant effects because in the large 
Nc limit the C/a(1) anomaly disappears and the rj' becomes also a pseudo-Goldstone boson [31] . 

Again we clarify that when fitting our theoretical formulas with the experimental data, all the masses 
of the pseudo-Goldstone bosons are taken from PDG [T], which we summarize in the Appendix [El for 
completeness. Only when discussing the Nc trajectories we start to use our prediction for the pseudo- 
Goldstone masses. The leading order mass parameters rri-^ = 2Bmq and rriK = B{mq + rus) in Eq.([T]), 
do not vary with Nc because B ~ 0{N^). This follows from the expression for the quark condensate in 
the chiral limit from Eq. ([T]), (0|^g-^|0) = —F'^B6^^, taking into account that both the quark condensate 
and F'^ are proportional to Nc- The bare masses are fixed in terms of the physical masses of the pion 
and kaon by employing the expressions given in Eqs. (IB.2P and (lB.4p . with the resulting values 

= 139.5^4:6 MeV , tTik = 519.6^7^5° MeV . (71) 
We summarize here the leading Nc scalings of the remaining parameters entering in our equations [40^ [48] 

A2~^, Mo~^^, {cd,Cjn,Cd,Cm,Gv} ^ y/N^ , 

l\c yNc 

{Ms, , Ms, , Mp , Mk' ,M^,M^, gsl} ~ 0{l) . (72) 
As commented above the U (3) large Nc relations between the singlet and octet scalar couplings 35] 
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are well fulfilled by our fitted values in Eq. ([55]) . while the mass relation Ms-^ = Msg is affected by a 30% 
variation, still within expectations for a leading large Nq prediction. 

The weak pion decay constant deserves a special attention because we can also calculate the 
subleading term in 1/Nc- The piont is that we have a 5-expansion for this important parameter in terms 
of F, the pseudoscalar weak decay constant in the chiral and large Nq limit appearing in the chiral 
Lagrangians, and that scales as y/Nc- The relation between F^ and F is given in Eq. ([02]) . Prom this 
expression the subleading 1/y/Nc contribution to the running with Nc of F.,^, is determined. This is, of 
course, a specific feature of employing U (3) xPT with its associated (5-expansion. 
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Figure 8: (Color online.) Masses for vr, r] and rj' as a function of Nc from 3 to 30 with one unit 
step. The different points are obtained by using the best fit given in Eq. (l55|) and the shadowed regions 
correspond to the error bands. Squares (magenta) are for rj' , diamonds (blue) for K, bursts (brown) for 
T] and circles (red) for vr. 



All the parameters in Eq. (f72l) . except asL, are present in chiral Lagrangians, namely in Eqs. ([I]), ([5]), 
(©, dl]), © and ([12]). It is important to notice that by construction, since U{3) xPT [331(391 BO] is 
a combined 1/Nc and chiral expansion, every coefficient multiplying a monomial of the fields does not 
contain any extra subleading piece in the I/Nq and chiral quark mass expansion. These extra terms 
would be part of higher order monomials in the (5-expansion. However, when fitting data, the actual 
numerical values in Eq. (|55p certainly reabsorb de facto higher order contributions in this expansion. 
This is also the case for the subtraction constant asL- 

We first consider the dependence on Nq of our results by taking the leading scaling with Nc of the 
parameters in Eq. ([72]) (but keeping the full Nc dependence of from Eq. ([C.2h . as follows from our 
calculation.) Later, when considering the properties of the a meson, we will also discuss variations in our 
results by considering subleading terms in the running with Nc for many of the parameters in Eq. ([72|) . 

Our full result for the Nc dependence of the masses of vr Eq. (IB. 21) . K Eq. ([B.4|) . ry and t]' Eq. ([T8]) . 
as well as the leading order mixing angle 9 Eq. ([B.7p . are shown in Figs. [8]and[9l respectively. The most 
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Figure 9: (Color online.) Leading order rj — rj' mixing angle 9 Eq. (IB .71) as a function of Nq from 3 to 30 
with one unit step. Similar to Fig. [8l the circle (red) points are from the best fit given in Eq.(|55p and 
the shadowed area corresponds to the error bands. 



streaking fact is the reduction by more than a factor of 2 of the mass of the rj with Nc- The mass of 
the 7]' also diminishes significantly. The K and vr masses vary little, specially the latter. For large Nc 
all the masses of the pseudo-Goldstone bosons n ,K ,rj and rj' go to zero in the chiral limit. However, for 
non-vanishing quark masses the rj' meson still gains a relatively large mass in the large Nc limit while 
the T] becomes similarly light as the pion [93]. This can be easily understood by looking at the leading 
order prediction in the large Nc limit for the pseudoscalar masses, 

ml = ml, 

m^/ = 2m|; — fnl , (74) 

where m.^ and m^/ stand for the masses of the rj and rj' mesons, corresponding to the notations of rrijj 
and m^' in Eqs. (|B.5p and (|B.6p with Mq — )• 0. Now, taking into account the values given in Eq. ([7T|) for 
and niK we then end with the following prediction for the leading order masses of rj and ry' in the 
large Nc limit 

rrirj = 139.5lt6 MeV , m^' = 721.5^^^^ MeV . (75) 

This leading order result already explains qualitatively the Nc behaviors for the masses of r] and rj' 
mesons shown in Fig. [8l 

For the leading order mixing angle 9, its Nc dependence is dominantly governed by the Ua_{^) anomaly 
mass Mq, as one can see from Eq. ()B.7p . It is also easy to demonstrate that in the large Nc limit, i.e. 
Mq = 0, the leading order mixing of r/ — r/' turns out to be the ideal mixing: 9 = —54.7°, as it should. 

It is worth stressing that previous works discussing the Nc behaviors of resonances, such as Refs.[T9l 
I21j . directly identify the t] meson with r/g, since it is based on SU{3) xPT. The rji is considered implicitly 
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through higher order counterterms, as any other heavy field in xPT. In addition, all of the pseudo- 
Goldstone masses are not changed in the variation of Nq [13 [201 [ZD [II]- However, as we have just 
shown, the change with Nc for the r] mass is very pronounced. 

Next we discuss our findings for the Nc dependence of the resonance properties channel by channel. 
As a clarifying remark let us mention that the following discussions are based mainly on the resonance 
poles appearing in the most relevant Riemann sheets for the energy region under discussion, as previously 
elaborated, unless a specific statement is given. 



Poles in the IJ = 00 channel 

We can easily track the pole trajectories for the cr, /o(980) and /o(1370) resonances while varying Nc- 
We first discuss the case of the a in Fig. [10] and latter we will turn our attention to the more massive 
/o(980) and /o(1370) resonances in Fig. [U] 
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Figure 10: (Color online.) Dependence of the a pole position, as a function of Nc from 3 to 30 with one 
unit step. In the left panel we show ^Js^ and in the right one Sa. The full result (circle in red) corresponds to 
our calculation without any approximation. The points labeled as vector reduced (triangle in blue) are obtained 
by keeping only the leading local terms from the vector resonance propagators. The empty squares (in green) 
correspond also to our full result but taking into account subleading terms in the running of Nc for the resonance 
couplings, according to the model of Eq. (|88l) . The shadowed areas enveloping the circles and the triangles show 
the error bands of our full result and vector reduced approximation. 



Fig. [To] has two panels. In the left one we show the running with Nc, starting with Nc = 3 in one 
unit steps, of the a pole position in the variable ^/s, similarly as done in Refs. [El [20l [95l [21], while 
in the right panel the same is shown employing the variable s, as in Refs. |961 I21| . Our full results 
correspond to the filled circles. The shadowed areas around every curve are generated by employing the 
same configurations of parameters that we exploit before in calculating the error bands from our fits in 
Eq. (155p and in Figs. [5][7] In this way we obtain the error bands for the shown Nc trajectories of the 
0" in Fig. [To] As one can see from this figure the resulting curves are quite stable even after taking into 
account the uncertainties of the inputs. From the left panel we observe that the width increases very fast 
with Nc, so that already for Nc = 7 it doubles as compared with its value at Nc = 3. In regards to its 
mass it first increases with Nc but for Nc above 7 it decreases. In the variable y/s the pole stays deep 
in the complex plane, as observed already in Refs. |19l [96l [2T1 124j . However, there is no point to keep 
interpreting the imaginary part of at the pole position as one half of the width when the mass is by 
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far much smaller. The correct interpretation, as pointed out in Refs. (961 121j . is obtained by considering 
the pole trajectory in the s complex plane: the mass square of the a resonance becomes real negative. 

In Refs. [2H the large Nc limit of the one- loop I AM is studied in the chiral limit and the following 
condition in terms of the 0{p'^) SU{3) chiral counterterms L2 and L3 [25] is obtained 

' > , So- approaches positive real s axis , 
25L2 + II-L3 < = , So- moves to 00 , (76) 
^ < , Sa approaches negative s axis . 



Making use of resonance saturation [35] 
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from the fitted values of Eq. (j55p . and taking My as Mp, we have 



- + 2MT ' ^^^^ 



25L2 + IIL3 = -11.2t°| X 10"^ < . (78) 

Thus, our curve with circles for Sa- in the right panel of Fig. [TOl bending towards negative large real 
values of s, is in agreement with the condition Eq. (|76p obtained in Ref. ^21j . 

Compared with the one-loop lAM result of Ref. [2Tj our curve bends faster to the left in the s-complex 
plane. It is interesting to point out that we can give rise to closer curves to those of Refs. [191 121] if 
we freeze out the full propagators of the vector resonances exchanged in the crossed channels and only 
take the leading local terms generated by them. This in fact corresponds to integrating out the vector 
resonance states in our theory and keeping only the leading contributions in the low energy sector, which 
are 0{p^). These terms are the ones incorporated in the one- loop xPT calculation used in NLO lAM 
Refs. [19] 121] . but not the higher order ones that arise by keeping the full vector resonance propagators 
that enter in our full calculation. The corresponding trajectory within this approximation is shown in 
Fig. dniby the triangles, labeled as vector reduced throughout. Notice that now the mass of the a keeps 
increasing with Nc, and the resulting trajectory is quite similar to that obtained in Ref. |19j . particularly 
for Nc ^ 10. Let us remark also that Refs. |19[ 195] . within their estimated uncertainties, also have cr 
pole trajectories in the y/s variable with decreasing mass (all of them has increasing width as our case.) 

One could also interpret the discrepancy between the circles and triangles in Fig. [10] arguing that 
higher order effects in the expansion employed could potentially have a large impact on the results. Let 
us first note that this criticism, from what we explicitly showed (the comparison between the triangles 
and circles in Fig. [TO]) . applies strictly to the lAM results. The higher orders included in our approach 
give rise to a large variation with respect to the results obtained when only their local reduction is kept, 
as done in the lAM. Second, from this one could of course infer the possibility that higher orders not 
considered in our approach, as well in any of other unitarization methods, could give rise to quite different 
a pole trajectories for large enough Nc- This is further discussed below, see also Ref. [95j . 

In contrast, to integrate out or keep the full contributions of the bare scalar resonance exchanges 
from £5 Eq. ([5|) does not change the Nc trajectory of a in a significant way. This also indicates that 
the bare scalar nonet has little influence on the a pole, as already pointed out in Ref. [10] . There it 
was stressed that the a, k and ao(980) originate independently of whether the bare scalar singlet and 
octet resonances are included in the formalism. This is also understandable by looking at how the 
resonances contribute to the LECs in xPT. For example, Li , L2 in SU{'i) are purely contributed 
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by the vector resonances, which also dominate L3 [l8j. For L4 , , Lq , Lg, they are dominated by the 
scalar resonances [l8]. Precisely in the vrvr — )• vrvr process, the factors accompanying L4 ,L5,Lq,Ls are 
always proportional to sm^, trn^ or m^, which are less important in the resonance region than the terms 
Li , L2 , that are accompanied by s^, st or t^. Indeed, in the chiral limit only the contributions to vrvr 
scattering proportional to Li, L2 and -L3 survive and with those LECs a and p were perfectly generated 
in Ref. [IHIM]. 

To conclude we can say that our results for the trajectory of the a resonance pole are in qualitative 
agreement with those from the one-loop JAM |191 [21} [M] and differences can be understood in terms of 
the higher orders included by using full resonance propagators for the crossed vector exchanges. 

Indeed, one-loop lAM is a particular case of our approach that results by expanding A^"^ = A^^^ (1 — 

one has 

T = Ar(2) (Ar(2) _ Ar(4) + Ar(2)^Ar(2)) . (79) 

In the previous expression the chiral order is indicated by the superscripts. For simplicity in the notation 
we have removed the subscript J and superscript / present in Eq. ([30]) . The derivations apply to any 
specific partial wave. Now, reducing Eq. ([30]) to 0{p^) and applying resonance saturation one has that 

iV(4) = + A(2)5iv(2) , 

T = r(2) rr(2) - r(^)j ""^ r^^) , (so) 

which is the one- loop lAM result. 

However, our results are quite different from those obtained in Refs. [211 [20] by employing the I AM 
with a two loop SU{2) xPT calculation for vrvr scattering, where for large Ac the a pole falls down to 
the positive real s-axis at around 1 GeV^. Indeed, as discussed in more detail below, we also obtain a 
pole in the large Nc limit at y/s = Ms^ — 1 GeV, but it comes from the bare singlet state which, for 
Nc = 3, is part of the /o(980) resonance. However, the a pole is not affected by its presence in a first 
approximation and this is not the reason for the different cr-pole trajectory obtained by us in comparison 
with Refs. [2ll[20]. 

In order to understand better this discrepancy with the two-loop lAM calculations [20] [2T] let us 
discuss the effects of several contributions particular to our approach. The first one corresponds to the 
influence of having U{2>) xPT. Compared with SU{2>) xPT, where 77 is identified with the octet r/g and r/i 
is only implicitly included at 0{p'^) through the LEG L-j , the present discussion incorporates the explicit 
contributions from the physical states of rj and r/' within U (3) with the 5 counting. As a result we 
take into account the drastic reduction of the r/ mass with Nc as discussed above, see Fig. [8| However, 
since the couplings of the a to ryry, 7777' and r/'r/' are weak in our calculation (but not necessarily in others 
|88il89| ). Tabled! one should then expect small effects on the cr-pole trajectory from this more exhaustive 
tTTataent of thi, and pseudoscalars in f/(3) xPtI To make this statement more quantitative it is 
interesting to find a way to make closer our approach to the 5'C/(3) one and then compare. This can 
be achieved by imposing the following conditions when varying Ac: (i) Freezing the mixing angle 
at leading order in Eq. ([13p to zero, together with the higher order mixing parameters in Eq. ([14p . (ii) 
The masses of the (then ?7g), vr , A are fixed to their physical values and do not vary with A^c- (hi) 



*This is also derived in Ref. [57] . 

^The weak couplings of the a to states with r) and r)' is kept when varying Nc as we have explicitly checked, even though 



the 77 mass decreases significantly with Nc , as shown in Fig. [Q] 
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The r]i mass is also fixed by using the leading order result from Eq.([T]) without mixing, which gives the 
value 1040 MeV. In this way, the r] meson in U{3) xPT resembles the rjs in SU{3) xPT and the rji only 
appears in the scattering amplitudes involving it , K ,r] mesons through the loops, which is suppressed 
by 1/Nc- This mimics the role of Lj in SU{3) xPT [19], although one needs to be careful about the 
Nq counting for Lj [971. The corresponding Sa Nq trajectory within this approximation, named as 
mimic SU(3), is so similar to our full results in Fig. [10] that we do not show it explicitly. Thus, for the 
fj case differences that arise by having used the U{3) xPT in our approach are not significant and are 
not certainly responsible for the differences with respect to the two-loop lAM results |21|, 120]. However, 
the situation could change dramatically for other resonances, e.g. for the /o(980) and /o(1370), and in 
general for any other resonance that had large couplings to states including the 77 and r]' mesons. We 
also want to point out that, by the same token, our study clarifies that previous results, e.g. those from 
Refs. \21 \ [20]. are stable under the explicit inclusion of the r/i singlet when increasing Nc- 

Another issue to be considered is the fact, pointed out repeatedly in Refs. [95j regarding the lAM, 
that the large Nc limit is a weakly interacting limit while the phenomenological success of unitarization 
methods, e.g. ours, typically rests on resumming infinite string of diagrams which are particularly 
strong. In our approach two types of resummations are considered simultaneously: (i) At the tree-level, 
by employing explicit resonance fields, so that infinite local terms are resummed by taking vector and 
scalar bare resonance exchanges, (ii) At the loop level, by resumming the RHC. The former resummation 
is leading in large Nc while the latter is subleading. Point (i) is the most important for studying the 
vector resonances, while point (ii) dominates for the scalar ones in real life [10 [ I95 [ IM]. 

Concerning the point (i), let us illustrate it by considering tree- level IJ =11 tttt scattering from the 
leading order Lagrangian in Eq.([T|), required by current algebra, and the vector resonance Lagrangians 
in Eqs.® and ([8]). It is well known the dominant role of the p resonance in IJ = 11 vrvr scattering at 
low and intermediate energies. These tree-level contributions are the ones that survive in the large Nc 
limit. The corresponding amplitude is 

-f 12M^^ + 24M25] + 6M2(4m2 - s - 1M1){2s + - Ami) l"g ^ ^ 4^2 | ' (^1) 

where the first term in the right hand side of the above equation is from the current algebra in Eq.([T|), 
the second term corresponds to the s-channel exchange of a p resonance and the last term is contributed 
by the p-exchanges in crossed channels, t and u. 

Let us now calculate the local terms at 0{p'^), T^'^^'^{s)2, and 0{p'^), T^^^°^{s)4 corresponding to 
the chiral expansion of Eq. (j8ip at low energy. We have for these terms: 
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An interesting task is to apply the lAM at NLO i5[[6l[71[2ll[2l] in order to reconstruct the full Tl'^'^°°{s) 
amplitude from the contributions in Eq. ([52]) . It results 
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Then, we see that only when 

0. = ^ ,84) 

is fulfilled, one can reproduce the exact input s = M^. The relation in Eq. (j84p is nothing but the so-called 
KSRF relation [98]. Within chiral resonance Lagrangians it was derived originally in Ref. [23] by studying 
vrvr scattering and later in Ref. [99j considering the pion vector form factor. We confirm this relation here 
from a different point of view. Due to the inclusion of the contributions from the crossed channels, the 
current version differs from the original one with Gy = Nevertheless if the contributions from the 

v2 

crossed channes are neglected, by dropping the last term in the right hand side of Eq. (|8ip . we can then 
recover the original KSRF relation Gy = ^ [49j 

Another look at Eq. (j83p can be obtained by employing the results of Ref. [24j which gives the p mass 
square in the chiral and large Nc limit from one-loop lAM at the position 

(85) 



4L3 



Taking into account Eq. (I77p (and keeping the vector contribution proportional to Gy which is much 
larger than the scalar one) one then has 

= ■ (86) 

Then only when Gy = Fj^/\/3 one can recover the exact position of s = My = M^, which is the same as 
Eq.dMD. 

This simple exercise shows that the success of the JAM in the vector channels [24] is intrinsically 
related with the KSRF relation and Vector- Meson-Dominance |100] . This success is not related with 
having resummed the RHC, as stated confusingly in Refs. [95], it occurs at the tree-level, and it is then 
expected to occur at any order in large Nc because loop corrections are even further suppressed with 
increasing A'^^. Notice that the benefits of the tree-level JAM resummation of Eq. (j83p is accomplished in 
our approach by the explicit exchange of bare vector and scalar resonances. Indeed, if Gy / F^/V3 the 
lAM resummation would not be so accurate, while including explicit resonance exchanges is not affected 
by tuning coupling or mass parameters. For an explicit example on this, corresponding to the exchange 
of a Higgs particle in the spontaneous breaking of electroweak symmetry, see Ref. [57] . 

Regarding the RHC resummation, point (ii) above, the situation is the opposite since loops in general, 
and unitarity loops in particular, become less relevant when Nq increases. In this respect, the situation 
is more worrisome than for the vector channels since when considering large values of Nc one is testing 
contributions not relevant in the physical Nc = 3 case. Within one-loop lAM it was shown in Ref. |2^ 
that this could imply inconsistencies. In Ref. [95] it is argued that if lAM describes data and resonances 
within a 10 to 20% errors, this means that the other contributions at Nc = 3 are not badly approximated. 
But since meson loops scale as 3/Nc, while tree-level inadequacies scale as 0(1), those 10 to 20% errors 
at Nc = 3 become 100% for values Nc ~ 30 and ~ 15, respectively. This criticism also applies to 
our approach because we rely in an expansion for the interacting kernel Nj{s), Eqs. (I30p and (136p . In 
this sense, at least, part of the differences between our results and two-loop lAM should correspond to 
different treatment of higher order terms in the chiral expansion, that at Nc = 3 are not so relevant 
but become more important for large Nc- Nevertheless, we think that the final reason why our results 
differ with respect to the two-loop lAM ones should be deeper, since already from the very start our 
curve for the full results on the right panel of Fig. [10] bends to the left while the analogous one for the 
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two-loop lAM results of Ref. 121] bends to the right for Nc > 3. There is also a big qualitative change by 
passing from one- to two-loop lAM results (which is not either properly understood yet in the literature 
and that is related to the issue we are discussing.) Our calculation certainly resembles much more to 
the one-loop case. On the one hand, we employ explicit resonance fields, with its associated tree-level 
resummation as just discussed. Note that the tree-level resummation for the lAM only generates the 
full resonance exchanges in the s-channel. E.g. taking again the toy model described above the lAM 
result, Eq. (f83|) . generates the p pole in the s-channel but no the cut due to p-exchanges in the crossed 
channels. The latter correspond to the log in Eq. (jSip . However, on the other hand, standard 
two-loop calculations [lOlJ. as those used in two-loop lAM ^211 [20] include more 0{p^) operators than 
those that can be reproduced by expanding the resonance propagators up to this order included in our 
study [51] (but let us recall that up to 0{p'^), and also up to 0{5), we have all of them.) The situation in 
regards to terms that become more important for the study of scalar resonances for large Nc, but that 
are not so important at Nc = 3, is then still far from being settled. 

Interestingly, in our approach one knows beforehand which is the resonance spectrum for Nc — oo. 
The reason is just because in the equation that we use for calculating the Tj-matrix, Eq. (pO]l . Nj{s)^^ 
scales as 0{Nc) while g{s) scales as 0{1). Then in the Nc — )• oo limit Tj{s) — )• Nj{s) and the latter 
only contains the poles included explicitly as bare resonances (vector and scalars) from the Lagrangians 
Eq. ([5]) and Eq. ([6|). In this sense, when reproducing data we are at the same time directly testing an 
Nc — 7" oo spectrum which is known beforehand. However, in the lAM due to the proliferation of higher 
order counterterms when passing from one- to two-loop order is not so clear which are the resonance 
poles that come up from the model for Nc — ?• oo, since they depend on the specific values of the chiral 
counterterms. 

Another issue of interest on which we want to elaborate further about its possible implications is 
the assumption of using the leading running with Nc for the parameters in Eq. (I72p . The possible 
impact of subleading terms in the running of Nc for chiral counterterms was already pointed out in 
Refs. [iHl [20l EU [24] ■ Refs. [El [20l [95] estimated this uncertainty by varying the xPT renormalization 
scale in which the counterterms are calculated. These values are then taken as the initial ones to engage 
the simple leading large Nc running. In our case, we include a subleading dependence on Nc on the 
coupling parameters. We think that for the bare masses this correction should be quite small because, 
as discussed below, the displacement of the bare resonance mass to its final pole position (at Nc = 3) 
is quite short. Since both limits are so close it is then a good approximation to keep them as 0(1) and 
not to elaborate further on subleading contributions for the bare masses. For the subtraction constant 
a^£, we do not consider any subleading term because at Nc = 3 it is already a common number, being 
the same for all the five channels with I, J = 00 as required by the large Nc U{3) symmetry. The 
parameters A2 ~ ^/Nc and Mq ~ l/^/Nc already vanish in the large Nc limit so that subleading terms 
are not so interesting because they are even further suppressed. Then, their contributions in the large 
Nc should be marginal compared with those that we keep. Coming back to the resonance couplings 
in Eq. (j55p we take into account additional QCD-inspired assumptions of high-energy behavior, such 
as unsubtracted dispersion relations for the pion electromagnetic form factor |23[ l99] and for the scalar 
strangeness changing scalar form factor |45[ I102j one has the relations, 



In the last equation the sum extends over the set of bare scalar resonances considered. In our approach we 
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only include one scalar octet of resonances, = 1. This is the so-called single resonance approximation 
which is based on the suppression for more massive resonances. The single resonance approximation 

was also employed in Ref. [M| for studying the p and a resonance pole trajectories within one-loop lAM 
in the chiral limit. Next we apply the large Nc relations of Eq. (j87l) to the physical case. It follows then 
the scaling with Nc'- 



Gv{Nc) = Gv{Nc = 3) 
Cd{Nc) = Cd{Nc = 3) 



CmiNc) = c^{Nc = 3) 



F^{Nc = 'i) 
F^Nc) 



FANc = 3) 



In the last equation we have assumed the same subleading dependence on Nc for Cm as for Cd- Note that 
Cm is poorly determined by the fit in Eq. (I55p and it is not excluded that Cm ~ Cd, as used e.g. in Ref. [45]. 
We also make use of the already discussed U{3) large Nc relations Cd = Cd/VS and Cm = Cm/V^ so that 
the scaling with Nc of all the couplings in Eq. ()55p is driven by F^^. Similar relationships were employed 
in Ref. [23] to deduce constraints between resonance saturation and unitarization by one-loop lAM. But 
note that here we now take advantage of the fact that we can calculate both the leading and subleading 
terms in the 1/Nc expansion of Eq. ()C.2p . within the (5-expansionl§ 
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We show by the empty squares in Fig.[TO]the resulting a trajectories taking into account the subleading 
terms in the running with Nc for the resonance couplings according to Eqs. (jSSp and (j89p . We observe 
that the curve still bends to the left, and qualitatively speaking is not very different to the circles 
obtained by considering the leading scaling with Nc shown in Eq. ()72p (we recall that we always include 
the subleading term in 1/Nc for F^, Eq. ([5^ .) It seems plausible then that subleading terms in the large 
Nc expansion are not either responsible for the difference between our results and two- loop lAM [2H |20] . 

The Nc trajectory for /o(980) is not reported in the lAM method due to the difficulty to track this 
pole [19]. In contrast, in our case its variation when moving Nc can be easily followed. The trajectories 
of the /o(980) and /o(1370) also pose a strong evidence that they are originated from the singlet and 
octet scalar states in the Lagrangian, respectively. However, one should also take into account that the 
contribution to the physical /o(980) due to a KK bound state O HOl [Ml IHU [92] disappears in the large 
Nc limit. We see that the difference in the trajectories obtained in the full result and mimic SU(3) cases 
differ significantly for the /o(980) and /o(1370) resonances. This is expected because these resonances 
couple strongly to states with the r] and tj' . For the /o(980) and /o(1370) resonances, since in our case 
their poles always fall down to the real axis at the masses of the singlet and octet bare resonances in the 
Lagrangian, in order, we only show the pole trajectories from the best fit in Fig. [TTJ For Nc — ?• oo the 
uncertainty is already given because their pole positions correspond to the bare masses of the singlet and 
octet scalar resonances, in order, given in Eq. (I55p together with their errors. 



However, note that in Eq. (|C.2|) we do not have all the counterterms of Cg2 because we use resonance saturation from 
the Lagrangians ([5]). Extra resonance operators would be needed which are beyond of our study [51j . 
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Figure 11: (Color online.) Prom top to bottom and left to right: Dependence on Nc from 3 to 30 with 
one unit step for the pole positions of the /o(908), /o(1370), Kq(1430) and ao(1450), respectively. The 
notation is the same as that in Fig. [101 For the last panel dedicated to the ao(1450) the Nq points for 
^ 7 are indicated explicitly for clarifying purposes. All the points with the same color correspond to 
the same framework (except the blue point in the ao(1450)-panel that corresponds to the mimic SU(3) 
case with Nc = 3.) 



Poles in the IJ = |0 channel 

The -ftrQ(1430) pole can be easily tracked in the complex plane, which we show in the third panel of 
Fig. [TTJ Like the /o(1370), this pole falls down on the real axis for Nc — t- oo, at s = MJ^, indicating that 
it originates from the bare octet scalar states with / = 1/2. At Nc = 3, its pole is somewhat sensitive to 
the different modifications of our results already discussed, while for large Nc they run towards the same 
point, as expected. From Tabled] we can see that the iro(1430) resonance couples to K-q' as strongly as 
to Ktt for Nc = 3, while its coupling to Kr] is suppressed. This behavior is also kept when varying Nc- 
For the k pole trajectories, depicted in Fig. W2\ the situation is more involved. After Nc = 4, we could 
not find its pole in the second Riemann sheet. However, switching to the fourth sheet, i.e. with (+, — , +), 
we can then track the k pole for any value of Nc- This is due to the more complicated cut structures in 
IJ = io case so that the trajectories of k are not as smooth as for the a in the IJ = 00 case. This is 
the reason why the k pole in the second Riemann sheet cannot be tracked because at Nc = 5 the pole 
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Figure 12: (Color online.) Nc dependence of the k pole position in the fourth Riemann sheet (+, — , +) 
from Nc = 3 to 30 with one unit step. The notation is the same as in Fig. [TOl 



moves to another Riemann sheet obtained by crossing the new cut in the complex plane, so that only a 
bump remains in the second Riemann sheet. Indeed, for Nq = 3 the k pole position in the fourth sheet 
is very similar to that in the second sheet for Nq = 3 (as discussed above.) The resulting trajectory in 
the fourth Riemann sheet is then shown by the circle points in Fig. [T2j We show the k pole trajectory 
with Nq both in the ^/s, left panel, and s variables, right panel. The same remark as for the a is in order 
also for the k. The correct interpretation is obtained by considering the pole trajectory in the s complex 
plane: the mass square of the k resonance becomes negative, as stressed in Ref. [96] for the similar a 
case. From Fig. [12] one observes that the pole trajectory bends to the left for the final as well as for 
the mimic SU(3) approximation. However, for the vector reduced approximation it tends to the right for 
Nc ^ 10 in the ^/s complex plane. In Fig. [12] the mimic SU(3) and vector reduced approximations are 
not shown, they were already shown in Fig. [TU] for the similar a case. 

Poles in the IJ = 10 channel 

The trajectories for the ao(1450) pole position are shown in the last panel of Fig. [TTl One can see that 
for this resonance the trajectories bend and make a knot for Nq < 7. In order to visualize this behavior 
different symbols for the pole positions up to Nc = 7 have been used in the figure. Independently from 
this peculiar behavior the ao(1450) pole position trajectories can be followed easily in a smooth way as 
Nc varies. As expected, for large Nc this pole moves to the real axis with zero width at s = MJ^, 
the same position as for the /o(1370) and ii'o(1430). It corresponds to the bare isovector members of 
the scalar octet Sg, Eq. ([7]). All SU{3) breaking effects in the masses of these heavier scalar resonances 
originate through pseudoscalar loops and disappear in the large Nc limit. 

For the ao(980) pole in the fourth sheet, we can track its trajectory up to Nc = 18 without any 
difficulty, which we display in Fig. (ll3p . After Nc = 19, only bumps appear in this sheet because there 
is a cut in the complex plane that connects the fourth Riemann sheet with another one, that would be 
obtained by crossing continuously this cut in the complex plane (not the real axis on which the standard 
Riemann sheets based on unitarity are defined), where the pole is finally located. Contrary to the k, case 
we cannot find a unique sheet where we can track the pole position from Nc = 3 onwards to high values 
of Nc- Although similar poles can be tracked if we switch to other sheets for Nc > 18, there is some 
discontinuity when passing from one trajectory to another in different Riemann sheets. The trajectory 
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Figure 13: (Color online.) Dependence on Nc for the ao(980) pole position from 3 to 18 with one unit 
step. The notation is the same as in Fig. [TOl 



shown in Fig. [13] corresponding to our fuh results clearly indicates that this resonance is dynamically 
generated (921 [9l [10] , since it moves further and further in the complex plane with a huge width increasing 
very fast with Nq- 

Poles in the IJ =11, 1 1 and 01 channels 

The Nc trajectories of the /9(770) and K*{S92) vector resonances are shown in the left and right panels, 
respectively, of Fig. [TJ] Both the p and K* move to the real axis with increasing Nc, as one can see in 
Fig. (|14p . indicating that their widths vanish in the large Nc limit. Indeed they vanish exactly as 1/Nc- 
In addition their masses move very little. Both facts are in agreement with its standard interpretation 
as qq resonances which implies an N^ scaling with Nc for the mass and 1/Nc for the width. 

Our results agree well with the previous conclusions on this respect [HI [20l [24]. Nevertheless, the 
Nc trajectories of the residues for the K*(892) in our full calculation show a clearly different feature, 
compared with the mimic SU(3) approximation. As one can see in Fig. [T5\ the residues of Ktt are 
rather similar between the two frameworks, while the residues of Krj' in the full result are obviously 
larger than the Krji in the mimic SU(3) approximation, which can be attributed to the rj — rj' mixing 
angle. In contrast to the almost flat residue of Krj (in fact Ktjq) in mimic SU(3) approximation, a more 
complicated structure appears in the full result. We verify this structure is caused by the Nc variation 
of the r] mass. As one can see the kink in the residue of Krj happens around Nc = 14, where precisely 
the Krj threshold becomes lighter than the i^*(892) mass. 
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Figure 14: (Color online.) Dependence on Nc of the />(770) and i^*(892) pole positions as a function of 
Nc from 3 to 30 with one unit step. The notation is the same as in Fig. [TUl 



The Nc trajectory for the </)(1020) pole position is completely analogous to those discussed for the 
p(770) and -fr*(892), but even simpler because it is an elastic channel and the width of the (/)(1020) is 
very small, Table [TJ This is why we do not dedicate a separate figure for this case. This type of trajectory 
is then in agreement with the standard qq interpretation for this resonance. Nevertheless, there is some 
movement in the (/)(1020) mass due to the small variation of the nearby KK threshold with Nq, see 
Fig.m 
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Figure 15: (Color online.) The Nq running of the various residues of the K*{892) resonance. The 
absolute value of the coupling to Kn is given by the circular points, while for the other states the ratio 
to the Ktt coupling is shown. Triangles apply to K-q {Kr]s) and squares to Kt]' (Krji). The left (right) 
panel corresponds to our full {mimic SU(3) ) results. 
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6 Conclusions 



In this work, we complete the first one-loop calculation of meson-meson scattering in U{3) xPT in the 
literature, that comprises simultaneously a chiral and an 1/Nc expansion, the so-called S counting. It 
is also the first work in which the Nc trajectories of resonance pole positions are studied taking into 
account the fact that the r]i becomes the ninth Goldstone boson in the large Nc limit. This has a large 
impact in the hadron spectrum because the mass of the r] pseudoscalar in large Nq decreases drastically 
and tend to be as light as the pion. 

In our one-loop U{3) xPT calculation for scattering amplitudes, instead of the local terms at 0(5'^), 
we include explicit resonance fields whose exchange generates these terms, as well as higher order ones. 
In order to compare with data, including the resonance region, we have unitarized the previous one- 
loop amplitudes projected in partial waves. For the unitarization we employ a non-perturbative scheme 
based on the N/D method, where the right-hand cut is resummed while the crossed-channel cuts are 
perturbatively treated as given by the input partial waves from U (3) xPT. We achieve a good reproduction 
of meson-meson scattering data from vrvr threshold up to energies between 1.2-1.6 GeV, depending on the 
particular partial wave. 

We have then studied the spectroscopy content of our solution by considering the poles and related 
residues of the different partial waves. Various resonance poles in the complex energy plane are then 
found, namely, for the a, /o(980), /o(1370), k, ^^5(1430), ao(980), ao(1450), p(770), K*{892) and (/)(1020) 
resonances. The pole positions agree remarkably well with the PDG values [Ij. The corresponding residues 
are calculated as well, which give us the coupling strengths of every resonance to the different channels. 
The couplings of the a resonances were studied in detail. We first discuss that the ratio of the couplings 
of the a to K^K~ and vr+vr". Our result for this quotient, 0.44^qq2, is in close agreement with the 
determination in Ref . ^86j , where its not-small size is interpreted as an indication in favor of the glueball 
content of the a. In our case, this resonance originates dynamically because of the vrvr interactions 
jnUnilJl]. We further compared our results for the couplings of the a to rjrj, ijij' and ij'r]' with previous 
results from Narison and collaborators. These authors obtain a coupling of the a to 7]r] that depends 
markedly on the value taken for the ratio of the couplings of the a to K^K~ and tt+tt", while ours is 
small. For the couplings to r/ry' and r/'ry' our calculation is in agreement with the upper bounds from 
Ref. [HH]. 

We pay special attention to the Nc dependences of the poles for these resonances as well as for their 
residues. Our approach fills a gap in the literature because previous studies in the literature did not take 
into account the Goldstone boson nature of the r/i field in the large Nc limit, due to the vanishing of the 
anomaly with large Nc- The trajectories for the pole positions are obtained by taking the leading 
scaling with Nc of the Ua{^) anomaly mass, A2, the subtraction constants agL and the bare resonance 
parameters. The scaling with large Nc of F-,^ also includes subleading terms determined from our one- 
loop calculation in C/(3) xPT. The Nc dependences of the pseudo-Goldstone masses have been taken 
into account for determining the resulting trajectories. It is specially remarkable the large reduction of 
the T] mass with Nc, that becomes similarly light as the pion for large Nc- We discussed in depth the 
particular case of the a resonance. We show that our results are in qualitative agreement with the ones 
of one-loop lAM. We obtain that Sq-, the pole position of the a resonance in the variable s, bends towards 
the real negative axis. This behavior is the expected one according to a relation obtained in Refs. |21ll24j . 
given the values of the large Nc leading 0{p^) counterterms L2 and L3 that follow from our calculation 
making use of resonance saturation. We also find that higher order terms arising within our approach by 
keeping the full resonance propagators in the crossed vector resonance exchanges give rise to a qualitative 
different behavior in the -y/s-complex plane for the a case for not too large Nc values {Nc ^ 7). This is 
easily understood because for increasing Nc the loop contributions from crossed diagrams tend to vanish 
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faster with Nq than those from the resonance exchange diagrams. In this way, the cancellation between 
both contributions found in Ref. [10] at Nc = 3 for the resonant scalar channels below 1 GeV, is spoilt 
for higher values of Nc and then there is more sensitivity to crossed channel dynamics. As a result, we 
find that the pole position for the a resonances has a mass that typically decreases with Nc, while this is 
not the case if the vector resonance exchanges are reduced to the leading local term contribution. Similar 
results are also obtained for the k resonance. 

The last point is one of the most problematic issues in determining the large Nc trajectory for the a 
resonance pole. As stressed in Refs. [95] one is sensitive in this case to terms that are not so significant 
at Nc = 3 so that a testing ground for the model is lacking. Particularly taking into account that the 
tendency for |so-| is to increase with Nc at around or above 1 GeV^, so that model dependence becomes 
more important. Indeed, our results are quite different from the two-loop lAM ones [211 120] where 
bends towards real positive values of s, falling down to the real axis at around 1 GeV^. Surely the most 
important reason for such disagreement has to do with the different way higher orders are treated within 
two-loop lAM and our approach. In addition, we have argued about other reasons. We have excluded 
that the explicit consideration from our side of the rji as a degree of freedom, by extending SU{2) or 
SU{3) xPT framework to U{3) xPT, is not responsible for such disagreement. The basic reason is that 
the a resonance within our approach couples weakly to states containing the r] and r]' and, then, is not 
so sensible to any improvement in the treatment related to those states. This is an interesting result 
by itself because then our calculation establishes that previous studies [13 [211 [20] , are stable under the 
explicit inclusion of the r]i singlet even when increasing Nc- However, this is not the case for all the 
other resonances studied that have large couplings to states with tj and/or r]' and that it is a necessary 
step forward to consider U{3) xPT. 

We also considered the possible presence of subleading terms in the large Nc running of the resonance 
parameters in the resonance Lagrangians by employing model results in the literature that are obtained by 
imposing QCD-inspired assumptions for the high energy behavior of some form factors. We then obtain 
that they scale as F^r, whose scaling we calculate within U{3) xPT at the one-loop level. Subleading terms 
in the resonance masses are not considered because they move very little from Nc = 3 (actual physical 
poles) to Nc = oo (bare masses). They are not either taken into account for parameters vanishing in the 
large Nc limit. We conclude tentatively that subleading terms are not responsible for the discrepancy of 
our results with those of two- loop lAM |2H 120], We also pointed out that within our approach we know 
beforehand the resonance spectrum in the Nc — >• oo limit, because it directly corresponds to the bare 
resonances introduced. In this way, when fitting data we are testing a known Nc — )■ oo set of resonances. 
In the lAM the knowledge of the Nc — )• oo resonance spectrum is more uncertain at the level of two-loops 
because of the proliferation of higher orders counterterms and the dependence of the precise spectrum 
with the explicit numerical values of the counterterms. 

Regarding the /o(980) we obtain that in the large Nc limit it tends to a zero width pole position 
corresponding to the bare singlet scalar resonance around 1 GeV. It is worth stressing that in other 
studies a zero width pole with mass of 1 GeV was already found for large Nc from the evolution of the 
a pole trajectory [211 EQ]. However, the pole at 1 GeV disappears in those studies for Nc = 3 and only 
the £7 remains [20]. In our case, both states remain. On the one hand, we have the bare scalar pole 
around 1 GeV that contributes to the /o(980) resonance pole in Nc = 3 while, on the other hand, the a 
resonance originates dynamically mainly from pion interactions. In addition, the strong KK interactions 
near threshold in S-wave gives rise to another strong contribution to the /o(980) as a KK bound state. 
This contribution disappears in the large Nc limit and only the bare singlet scalar state pole contributes 
then. It is an interesting exercise for future work to check that our bare singlet scalar pole at around 
1 GeV is enough to guarantee local-duality [103j . 
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The ao(980) resonance disappears deep in the complex plane for large Nq, by increasing its width and 
mass, which is always positive and large. This behavior corresponds to a dynamically generated resonance. 
We also discussed the Nc trajectories for the /o(1370), ao(1450) and K*{1430). Asymptotically for large 
Nc, they tend to the zero width pole position of the bare octet of scalar resonances included at the 
tree- level around 1.4 GeV. 

Finally, the vector resonances p{770), K*{892) and (/>(1020) are reproduced with properties in good 
agreement with the PDG [I]. They have an Nc pole position trajectory in good agreement with the 
expectations for a qq state, with a quenched mass running as A'^^ and the width decreasing as 1/Nc- The 
decreasing r] mass with Nc makes that the K* (892) resonance becomes heavier than the K-q channel for 
Nc around 15. The crossover between the mass of the resonance and the Kr] threshold manifests in a 
kink in the coupling to Kr]. 
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Appendix 



A Convention of the loop functions 

The one-loop functions that appear in our ?7(3) xPT calculation are calculated in dimensional regular- 
ization within the MS — 1 renormalization scheme ^5]. They are defined as 



iTT^ J q^ — 
2 

-m in — , 
(2vrM)4-^ f^^^ 1 1 

2 — — p)2 — 77j2 



= 1 - log^ + x+log— t +x_log— , (A.l) 

/x^ x+ x_ 

where /i is the renormalization scale, s = and x± was defined in Eq. (l33p . In case of equal masses, the 
two-point function reduces to 

Bo{s, rr?) = l- log ^ + a{s) log "^'l ~ \ , (A.2) 
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with a{s) defined in Eq. ([35]) . 



B Renormalization of the wave functions and masses 

Prom the calculation of the tt, K,f] and fj' self-energies Fig. [21 with the latter fields defined in Eq. (fTH]) . 
we have: 



1 



2 2 
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8cm(cm - Cd)m|^(2m^ + m2) 8cm(cm - m^(m^ - m2) 
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3F2M| 



(B.l) 



(B.2) 



(B.3) 



The leading order masses of ry,/?', i.e. and m^/ defined in Eq. ()14p . and the 7] — rj' mixing 
from Eq.([T|) are found to be 



(B.4) 
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(B.5) 
(B.6) 
(B.7) 



with A2 = — ml . 
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The higher order mixing parameters of ry — ry' are defined in Eq. ()14p . They read 
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C Renormalization of pion decay constant 

The definition of pseudoscalar weak decay constant is 

(0|A^|/(p))=^F^.p^5'^^ (C.l) 

where the axial- vector current is A'^ = 97^75^(7. 

Throughout this work we have expressed F in the chiral hmit in terms of the physical pion decay 
constant The corresponding expression coincides with the one in SU{3) xPT [251 [38], once resonance 
saturation of the 0{p^) xPT counterterms is assumed [l8]. It reads 



FU + 



1 



167r2F2 ^ 

AcdCm{ml + 2m1) 8cd Cm {m\ - ml) 



+ 



3F2M| 



(C.2) 



D Scattering amplitudes 



By using crossing symmetry and isospin symmetry, as we discussed in Section [H all the meson-meson 
scattering amplitudes with well defined isospin and angular momentum in t/(3) xPT can be reduced to 
the calculation of 16 independent processes. 

Due to the much lengthy expressions of the one-loop U{3) xPT scattering amplitudes, we only provide 
here the analytical expressions for the tree level results at leading order, denoted with the superscript (2). 
For the remaining parts, comprising the ones from the loops and resonance exchanges, one can download 
the Mathematica file from http:/ /www.um.es/oller/u3FullAmpl6.nb^ that can also be provided by the 
authors under request. 
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E Numerical inputs 

In the numerical discussion, we take the average values for the masses of the charged and neutral pions 
and kaons. The values taken are summarized below in units of MeV: 

= 137.3 , rriK = 495.6 , = 547.9 , m^/ = 957.7 , 

= 92 A , n = 770.0 , = 1026.0, = 783.0 . (E.17) 

The value for is adjusted so as to reproduce the pole position of the (/)(1020) resonance in Table [T] 
according to the PDG [T]. 
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